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ABSTRACT 

Building on results from the Magnetism in Massive Stars (MiMeS) project, this 
paper shows how a two-parameter classification of massive-star magnetospheres in 
terms of the magnetic wind confinement (which sets the Alfven radius Ra) and stellar 
rotation (which sets the Kepler co-rotation radius Rk) provides a useful organisa- 
tion of both observational signatures and theoretical predictions. We compile the first 
comprehensive study of inferred and observed values for relevant stellar and magnetic 
parameters of 64 confirmed magnetic OB stars with T g > 16 kK. Using these param- 
eters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot 
of i?K vs. i?A- This diagram can be subdivided into regimes of centrifugal magneto- 
spheres (CM), with Ra > Rk, vs. dynamical magnetospheres (DM), with Rk > Ra- 
We show how key observational diagnostics, like the presence and characteristics of Ha 
emission, depend on a star's position within the diagram, as well as other parameters, 
especially the expected wind mass-loss rates. In particular, we identify two distinct 
populations of magnetic stars with Ha emission: namely, slowly rotating O-type stars 
with narrow emission consistent with a DM, and more rapidly rotating B-type stars 
with broader emission associated with a CM. For O-type stars, the high mass-loss 
rates arc sufficient to accumulate enough material for line emission even within the 
relatively short free-fall timescale associated with a DM: this high mass-loss rate also 
leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the 
longer confinement of a CM is required to accumulate sufficient emitting material 
from their relatively weak winds, which also lead to much longer spindown timescales. 
Finally, we discuss how other observational diagnostics, e.g. variability of UV wind 
lines or X-ray emission, relate to the inferred magnetic properties of these stars, and 
summarise prospects for future developments in our understanding of massive-star 
magnetospheres . 

Key words: stars: magnetic fields - stars: early-type -circumstellar matter - stars: 
mass-loss- stars: rotation- stars: fundamental parameters - stars: emission-line, Be - 
ultraviolet: stars - X-rays: stars. 
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1 INTRODUCTION 

Building on pioneering detections of strong (kG) fields in the 
chemically peculiar Ap and Bp stars (e.g. Babcock 1947; 
Borra & Landstreet 1980), new generations of spectropo- 
larimeters have directly revealed large-scale, organised (of- 
ten predominantly dipolar) magnetic fields ranging in dipo- 
lar strength 1 from order of 0.1 to 10 kG in several dozen OB 
stars (e.g. Donati et al. 2002, 2006a; Hubrig et al. 2006; Pe- 
tit et al. 2008; Grunhut et al. 2009; Martins et al. 2010). In 
recent years, an observational consortium known as MiMeS 
(for Magnetism in Massive Stars) has been carrying out sur- 
veys to detect new magnetic OB stars, while also monitoring 
known magnetic OB stars with high resolution spectroscopy 
and polarimetry (Wade et al. 2011a). Concurrently, theoreti- 
cal models (Townsend et al. 2005, 2007) and magnetohydro- 
dynamical (MHD) simulations (ud-Doula & Owocki 2002; 
ud-Doula et al. 2008, 2009) have explored the dynamical in- 
teraction of such fields with stellar rotation and mass loss, 
showing for example how suitably strong fields can channel 
the radiatively driven stellar wind outflow into a circumstel- 
lar magneto sphere. This paper aims now to provide an initial 
classification of the observed magnetospheric properties for 
a broad sample of magnetic massive stars. 

The idea of a magnetosphere has been exploited to ex- 
plain particular properties of some massive stars, for ex- 
ample the photometric light curve and Ha variations of 
the He-strong star a OriE (Landstreet & Borra 1978), the 
UV resonance line variations of magnetic Bp stars (Shore 
& Brown 1990), the X-ray properties of the O-type star 
9 1 OriC (Gagne et al. 2005), and the radio emission of Ap- 
Bp stars that correlates with the field strength (Linsky et al. 
1992). 

For a few specific stars, previous work has already 
shown some promising agreement between theoretical pre- 
dictions and key observational characteristics. For example, 
the luminosity, hardness, and rotational modulation of X- 
rays observed in the O-type star S 1 OriC all match well 
the X-rays computed in MHD simulations of magnetically 
confined wind shocks, which result from the collision of the 
wind from opposite footpoints of closed magnetic loops in 
its ~ lkG dipole field (Gagne et al. 2005). In the B2p star 
crOriE, the combination of its very strong (~ 10 kG) field 
and moderately fast (1.2-day period) rotation leads to for- 
mation of a centrifugally supported magnetosphere with ob- 
served, rotationally modulated Balmer line emission reason- 
ably well explained within the Rigidly Rotating Magneto- 
sphere model (RRM; Townsend & Owocki 2005; Townsend 
et al. 2005). Most recently, Sundqvist et al. (2012) showed 
that, even in the very slowly rotating (537-day period) O- 
type star HD 191612, the magnetic confinement and tran- 
sient, dynamical suspension of its strong wind mass loss 
leads to sufficient density to likewise match the observed 
rotationally modulated Balmer line emission. 

Building on these results, along with those from the 
MiMeS observational survey, this paper compiles an exhaus- 
tive list of confirmed magnetic, hot OB stars, along with 
their physical, rotational and magnetic properties (§2). As 



a basis for organising this compilation according to mod- 
elling predictions, we follow (§3) the two-parameter theo- 
retical study of ud-Doula et al. (2008), which characterised 
MHD simulation results according to the strength of wind 
magnetic confinement (77*) and fraction of stellar rotation 
to orbital speed at the stellar equatorial radius (W). These 
dimensionless parameters uniquely define associated char- 
acteristic radii, namely the Alfven radius Ra and Kepler 
co-rotation radius _Rk- 

We show in particular (§4) that an associated log-log 
plot of known magnetic stars in the Ra-vs.-Rk (or equiv- 
alently n^-vs.-W) plane, the magnetic confinement-rotation 
diagram, provides a particularly useful initial classification 
for interpreting the Ha properties of their associated magne- 
tospheres. Furthermore, we also explore the UV and X-ray 
characteristics as potential additional proxies of magneto- 
spheres (§5). We briefly review our main findings and con- 
clusions in §6. 



2 EXHAUSTIVE LIST OF MAGNETIC 
O-TYPE AND EARLY B-TYPE STARS 

A central goal of this paper is to compile a comprehensive list 
of OB stars for which magnetic fields have been convincingly 
detected via the Zeeman effect, so that their magnetospheres 
can be classified. 

The work here is done within the context of the MiMeS 
project (Wade et al. 2011a), which aims to expand the 
population of known magnetic stars, confirm the detection 
of poorly studied magnetic OB stars, and provide a mod- 
ern determination of their magnetic field characteristics. 
These goals are being achieved through Large Program ob- 
serving allocations at the Canada-France-Hawaii Telescope 
(CFHT), the Telescope Bernard- Lyot (TBL) and the ESO 
3.6m Telescope to collect high resolution, high signal-to- 
noise ratio spectropolarimetric observations of massive stars 
(see Wade et al. 2011a; Oksala et al. 2012; Alecian et al. 
2011, respectively). 

Table 1 lists our derived sample of 64 magnetic stars, 
ordered by spectral type and temperature. Column (1) gives 
the numerical identification (ID) we use in the figures. Col- 
umn (2) gives the HD number or a SIMBAD-friendly des- 
ignation. A dagger indicates that a note for that particular 
star is available in Appendix A. Columns (3) and (4) give 
a commonly used designation and the spectral type, respec- 
tively. Column (5) indicates if the star is a known single- or 
double-line spectroscopic binary (SB1/2), slowly pulsating 
B-type star (SPB), j3 Cep-type pulsator (/3Cep), or a Her- 
big Be star (HeBe). Table 2 compiles, for each star, the list 
of references where information can be found or how it is 
derived from MiMeS observations or other archival data. 



1 In the following all field strengths will be given as dipolar, un- 
less explicitly noted otherwise. 
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Tabic 1: List of magnetic massive OB stars and their physical, rotational and 
magnetic properties. 



ID 


Star 


Alt. name 


Spec, type 


Remark 


T e g 


logs 


log(L*/L Q ) 


R* 


M» 


P 


vsini 


Bp 












kK 


cgs 




Rq 


Mq 


d 


km s _1 


kG 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


HD 148937 




06 f?p 




41±2 S 


4.0 ±0.1 


5.8±0.1 


15 


60 


7.0323 


< 45 


1.0 


2 


CPD-28 2561 




06.5f?p 




35±2 S 


4.0 ±0.2 


5.5 ±0.2 


14 


43 


70 




>1.7 


3 


HD 37022 t 


6 1 OriC 


07Vp 


SB1 


39 ± I s 


4.1 ±0.1 


5.3±0.1 


9.9 


45 


15.424 


24 


1.1 


4 


HD 191612 1 




06f?p-08fp 


SB2 


35 ± I s 


3.5 ±0.1 


5.4 ±0.2 


14 


30 


537.2 


< 60 


2.5 


5 


NGC 1624-2 




06.5f?cp-08f?cp 




35±2 S 


4.0 ±0.2 


5.1 ±0.2 


9.7 


34 


158.0 


< 3 


>20 m 


6 


HD 47129 1 


Plaskett's star 


07.5111 


SB2 


33±2 S 


4.1±0.1 


5.09 ±0.04 


10 


56 




305 


>2.8 


7 


HD108 




08f?p 




35±2 S 


3.5±0.2 


5.7±0.1 


19 


43 


18000 


< 50 


>0.50 


8 


ALS 15218 t 


Trl6-22 


08.5 V 




34 ±2 


4.0 ±0.2 


5.0 ±0.1 


9.0 


28 




25 


>1.5 


9 


HD 57682 




09 V 




34 ± 1 s 


4.0 ±0.2 


4.8 ±0.2 


7.0 


17 


63.571 


15 


1.7 


10 


HD 37742 t 


C Ori Aa 


09.51b 


SB2 


29 ± I s 


3.2 ±0.1 


5.6±0.1 


25 


40 


7.0 


110 


0.060 


11 


HD 149438 


T SCO 


B0.2V 




32 ± I s 


4.0±0.1 


4.5 ±0.1 


5.6 


11 


41.033 


5 


0.20 m 


12 


HD 37061 1 


NUOri 


B0.5V 


SB2 


31.0 ±0.5* 


4.2 ±0.1 


4.4 ±0.1 


5.7 


19 




225 


0.65 


13 


HD 63425 




B0.5V 




29 ± I s 


4.0±0.1 


4.5 ±0.4 


6.8 


17 




< 10 


0.46 


14 


HD 66665 




B0.5V 




28 ± I s 


3.9±0.1 


4.2 ±0.5 


5.5 


9.0 


21 


< 10 


0.67 


15 


HD 46328 


£ x CMa 


B1III 


P Cep 


27±2 S 


3.5 ±0.2 


4.6 ±0.1 


8.6 


9.0 


4.26 


< 15 


>1.5 


16 


ALS 8988 


NGC 2244 OI 201 


Bl 


HeBe 


27±l s 


4.18 ±0.06 


4.05 ±0.08 


4.7 


12 




23 


>1.5 


17 


HD 47777 


NGC 2264 83 


B1III 


HeBc 


27±2 S 


4.0 ±0.2 


4.1±0.1 


5.0 


9.0 




65 


>2.1 


18 


HD 205021 1 


/3Cep 


B1IV 


SB2,/3Cep 


26 ± I s 


3.7±0.1 


4.22 ±0.08 


6.5 


12 


12.00092 


27 


0.36 


19 


ALS 15211 1 


Trl6-13 


Bl V 




26 ±2 


4.0 ±0.2 


4.0 ±0.1 


4.9 


9.0 






>1.4 


20 


HD 122451 1 


/3Cen 


Bl 


SB2,/3Ccp 


25 ±2 


3.5 ±0.4 


4.4 ± 0.2 p 


8.7 


8.8 




75 


>0.25 


21 


HD 127381 


a Lup 


B1/B2 V 




23 ± I s 


4.0±0.1 


3.76 ±0.06 


4.8 


9.0 


3.0197 


68 


0.50 


22 


ALS 3694 


NGC 6193 17 


Bl 




20 ±3 


4.0 ±0.4 


3.7±0.3 p 


5.6 


11 




83 


>6.0 


23 


HD 163472 


V2052Oph 


B1/B2V 


/3Ccp 


25 ± I s 


4.2 ±0.1 


3.8±0.1 


4.1 


10 


3.638833 


68 


0.40 


24 


HD 96446 1 


V430 Car 


BlIVp/B2Vp 




21.6±0.8 S 


4.0±0.1 


3.6 ±0.2 


4.5 


8.0 


5.73 


3 


6.5 


25 


HD 66765 




B1/B2V 




20 ±2 


3.9 ±0.2 


3.6 ± 0.2 c 


5.3 


7.5 


1.61 


100 


>2.1 


26 


HD 64740 


HR3089 


B1.5Vp 




24 ± 1 s 


4.0±0.1 


4.1 ±0.3 


6.3 


11 


1.33026 


160 


16 


27 


ALS 15956 


Col 228 30 


B1.5V 




23 ±3 


3.6 ±0.2 


4.3 ± 0.2 c 


9.1 


11 






>1.5 


28 


ALS 9522 


NGC 6611 W601 


B1.5Ve 


HeBe 


22±2 S 


3.8 ±0.3 


4.0 ±0.1 


6.4 


10 




190 


>4.0 


29 


HD 36982 


LP Ori 


B1.5Vp 




22 ±2 


4.0 ±0.2 


3.1 ±0.2* 


2.5 


2.2 




80 


0.91 


30 


HD 37017 1 


V 1046 Ori 


B1.5-2.5IV-Vp 


SB2 


21 ±2 


4.1 ±0.2 


3.4 ±0.2 


3.9 


7.2 


0.90119 


90 


>6.0 


31 


HD 37479 


o-OriE 


B2 Vp 




23 ±2 


4.0 ±0.2 


3.6±0.2P 


3.9 


5.0 


1.1908 


170 


9.6 m 


32 


HD 149277 f 




B2IV/V 


SB2 


22 ±3 


4.0 ±0.4 


4.0±0.3 p 


7.0 


17 




15 


>4.7 


33 


HD 184927 


V1671Cyg 


B2 Vp 




22 ± I s 


3.9±0.1 


3.6 ±0.2 


4.3 


5.5 


9.530 


14 


3.9 m 


34 


HD 37776 1 


V 901 Ori 


B2Vp 




22 ±1 


4.0±0.1 


3.5 ±0.1 


3.8 


5.5 


1.538756 


95 


15 m 


35 


HD 136504 1 


e Lup 


B2 IV-V 


SB2,/3Cep 


22 ±2 


4.0 ±0.2 


3.8±0.2 C 


5.3 


8.6 




42 


>0.60 


36 


HD 156424 




B2V 




22 ±3 


4.0 ±0.3 


3.7±0.4P 


4.8 


8.5 




15 


>0.65 


37 


HD 156324 




B2V 




22 ±3 


4.0 ±0.3 


3.7±0.4P 


5.1 


9.4 




60 


>1.8 


38 


HD 121743 


(j> Cen 


B2IV 


/3Ccp 


22 ±3 


4.0 ±0.3 


3.7±0.2P 


4.7 


8.0 




80 


>0.53 


39 


HD 3360 


CCas 


B2IV 


SPB 


20.4 ± 0.9 s 


3.8±0.1 


3.7±0.2 


5.9 


8.3 


5.37045 


17 


>0.34 


40 


HD 186205 + 




B2Vp 




20 ±3 


4.0 ±0.2 


3.5 ± 0.2 c 


4.9 


7.4 




5 


>1.7 


41 


HD 67621 




B2IV 




19 ±3 


4.0 ±0.3 


3.3±0.2P 


4.1 


6.2 


3.59 


20 


>0.90 


42 


HD 200775 t 


V3780 Cep 


B2Ve 


SB2, HeBe 


18 ±2 


3.4±0.2 


4.0 ±0.3 


10 


10 


4.328 


26 


1.0 


43 


HD 35912 


HR1820 


B2V 




18 ±1 


4.0 ±0.1 


3.3±0.3 p 


4.4 


7.2 


0.89786 


< 12 


>6.0 
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Table 1: Continued 



ID Star Alt. name Spec, type Remark T e g logg log(L*/Z/0) R* M* P vsmi B p 

kK cgs R Q Mq d km s" 1 kG 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(?) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


44 


HD 66522 




B2 III 




18 ±2 


4.0 ±0.4 


3.3±0.2P 


4.6 


7.6 




< 10 


0.90 


45 


HD 182180 


HR7355 


B2Vn 




17±l s 


4.2±0.2 


3.0 ±0.1 


3.7 


6.0 


0.5214404 


310 


11 


46 


HD 55522 


HR2718 


B2 IV/V 




17.4 ± 0.4 s 


4.2 ±0.1 


3.0 ±0.1 


3.3 


5.5 


2.729 


70 


>2.6 


47 


HD 142184 


HR5907 


B2V 




17±l s 


4.3±0.1 


2.8±0.1 


3.1 


5.5 


0.50828 


290 


10 


48 


HD 58260 1 




B3Vp 




20 ±2 


3.5±0.2 


4.1±0.2 C 


9.5 


9.5 




< 12 


>7.0 


49 


HD 36485 f 


5 0riC 


B3Vp 


SB2 


20 ±2 


4.0 ±0.1 


3.5 ± 0.1 c 


4.5 


7.1 


1.47775 


32 


10 


50 


HD 208057 1 


16 Peg 


B3V 


SPB 


19 ±3 


3.9 ±0.2 


3.6 ± 0.2 c 


5.5 


7.1 


1.441 


104 


>0.50 


51 


HD 306795 


NGC 3766 MG170 


B3V 




18 ±2 


3.9 ±0.2 


3.2±0.3 p 


4.1 


4.3 




65 


>5.0 


52 


HD 25558 1 


40Tau 


B3V 


SPB 


17±2 


4.0 ±0.2 


3.1 ±0.2* 


3.9 


5.5 




28 


>0.15 


53 


HD 35298 f 




B3Vw 




16 ±2 


3.8 ±0.2 


3.2 ± 0.2 c 


5.5 


5.6 


1.85336 


260 


>9.0 


54 


HD 130807 


oLup 


B5 




18 ±2 


4.1 ±0.1 


3.1±0.1 c 


3.5 


5.7 




25 


>2.0 


55 


HD 142990 t 


V913Sco 


B5V 




17±2 


4.2±0.2 


2.9±0.2P 


3.1 


5.7 


0.97907 


125 


>7.5 


56 


HD 37058 1 


V359 Ori 


B3VpC 




17±2 


3.8 ±0.2 


3.5 ± 0.2 c 


5.6 


6.6 


14.61 


25 


>3.0 


57 


HD 35502 1 




B5V 


SB2 


16 ±2 


3.8 ±0.2 


3.3 ± 0.2 c 


5.7 


5.7 


0.85 


80 


>6.8 


58 


HD 176582 


HR7185 


B5IV 




16 ± I s 


4.0±0.1 


2.9±0.1 


3.6 


4.7 


1.581984 


105 


7.0 


59 


HD 189775 


HR7651 


B5V 




16 ±2 


3.8 ±0.2 


3.2 ± 0.2 c 


5.3 


5.5 


2.6048 


85 


>4.5 


60 


HD 61556 1 


HR2949 


B5V 




15 ±2 


4.0 ±0.3 


2.6 ±0.11" 


2.8 


2.9 


1.9093 


70 


4.0 m 


61 


HD 175362 t 


Wolff's star 


B5V 




15 ±3 


3.7 ±0.2 


3.2 ± 0.1 c 


5.8 


5.3 


3.6738 


35 


>21 m 


62 


HD 105382 t 


HR4618 


B6III 




17±2 


4.0 ±0.2 


3.0±0.2P 


3.6 


4.8 


1.285 


90 


2.3 


63 


HD 125823 


aCen 


B7IIIp 




19 ±2 


4.0 ±0.2 


3.2±0.1 p 


3.6 


4.7 


8.812 


15 


>1.3 


64 


HD 36526 


V1099Ori 


B8Vp 




16 ±3 


4.0 ±0.3 


2.5 ±0.3? 


2.4 


2.0 


1.5405 




>10 



T Notes in Appendix. ~ — ~ 

(5) Single- or double-lined spectroscopic binary (SB 1-2), slowly pulsating B-type star (SPB), Cep-type pulsator (/3Cep), Herbig Be star (HeBe). 
s Parameters determined from modern spectral modelling. 

p Luminosity derived from our photometric calculations with bolometric correction. 
c Luminosity derived from SED fitting with Chorizos. 



Higher multipole components. 



Classification of Massive Star Magnetospheres 



Table 2: List of references for propeties retrieved from the literature with su- 
perscript letters indicating the type of parameter. Superscript letters in column 
(2) indicate information we inferred from MiMeS observations or other archival 
data, as indicated in the text. 



ID 


Star 


Ref. 


ID 


Star 


Ref. 


(1) 


(2) 


(3) 


(1) 


(2) 


(3) 


1 


HD 148937 


Wade et al. (2012b) s ™" 


32 


HD 149277™ 


Bagnulo ct al. (2006) 6 






A. Fullterton (priv. com.)" 






Landstreet et al. (2007) s 






Naze ct al. (2012b)* 


33 


HD 184927 


Wade et al. (1997) r 


2 


CPD-28 2561 


Barba et al. (MiMeS in prep) sr6a 






Yakunin et al. (MiMeS in prep) s6 "" 


3 


HD 37022 


Simon-Dfaz et al. (2006) a 


34 


HD 37776 


Landstreet et al. (2007) s 






Stahl et al. (2008)™ 






Kochukhov et al. (2011) 6 






Wade et al. (2006) 6 






Mikulasek et al. (2011) r 






Walborn & Nichols (1994)" 






Shultz et al. (MiMeS in prep)" 






Stelzer et al. (2005)* 






Shore & Brown (1990)" 


4 


HD 191612 


Wade et al. (2011b) a ™" 


35 


HD 136504" 


Uytterhoeven et al. (2005) sr 






A. Fullerton (priv. com.) 11 






Shultz et al. (2012) 6 






Naze et al. (2007)* 






Hubrig et al. (2009) 6 


5 


NGC 1624-2 


Wade et al. (2012a) a ™ a * 


36 


HD 156424 


E. Alecian (MiMeS in prep) s ™" 


6 


HD 47129 


Linder et al. (2008) srau 


37 


HD 156324 


E. Alecian (MiMeS in prep)*"' 6 " 






Grunhut et al. (2012b) 6 


38 


HD 121743 


Wolff (1990) s 


7 


HD108 


Martins et al. (2010) s ™ 






E. Alecian (MiMeS in prep)™" 






Marcolino et al. (2012)"" 






Grillo et al. (1992)* 






Linder et al. (2006)* 


39 


HD 3360 


Neiner et al. (2003) sr6tm 


8 


ALS 15218 


Gagne et al. (2011) s 






Oskinova et al. (2011)* 






Naze et al. (2012a)™" 


10 


HD 186205 


Zboril & North (2000) r 






Naze et al. (2011)* 






J. Grunhut (MiMeS priv. com.) s6 " 


9 


HD 57682 


Grunhut ct al. (2009) s6a "* 


41 


HD 67621 


Alecian et al. (MiMeS in prep) 3 ™" 






Grunhut et al. (2012c) r 


42 


HD 200775 


Alecian et al. (2008a) s ™ 


10 


HD 37742 


Bourct ct al. (2008) s ™" 






Hamaguchi et al. (2005)* 






Kaper et al. (1996)" 


43 


HD 35912 


Simon-Diaz (2010) s " 






Raassen et al. (2008)* 






Bychkov et al. (2005) r6 


11 


HD 149438 " 


Simon-Dfaz et al. (2006)" 


44 


HD 66522 


Zboril ct al. (1997) s 






Donati ct al. (2006b)™" 






Leone et al. (1997) s 






Mewe et al. (2003)* 






E. Alecian (MiMeS in prep) r6a 


12 


HD 37061 " 


Simon-Dfaz et al. (2011) a ™ 


15 


HD 182180 


Rivinius ct al. (2012) 3 ™" 






Petit et al. (2008) 6 


16 


HD 55522 " 


Briquet et al. {20QA) Br 






Stelzer et al. (2005)* 






Briquet et al. (2007) 6 


13 


HD 63425 


Petit et al. (2011) a ™ a " 


47 


HD 142184 


Grunhut ct al. (2012a) srba 


14 


HD 66665 


Petit et al. (2011) a ™ a " 






Oskinova et al. (2011)* 


15 


HD 46328 


Fourtune-Ravard ct al. (2011) s ™"" 


48 


HD 58260 


Bohlender (1989) s 






Oskinova et al. (2011)* 






Cidale et al. (2007) si> 


16 


ALS 8988 


Alecian et al. (2008b) s ™ 






Pedersen (1979)™ 






Wang et al. (2008)* 






Shore & Brown (1990)" 


17 


HD 47777 


E. Alecian (priv. com.)'™ 


19 


HD 36485 


Leone et al. (2010) 8 ™" 






Naze (2009)* 






Shore & Brown (1990)" 


18 


HD 205021 


Donati et al. (2001) s ™" 


50 


HD 208057 


Chauvillc et al. (2001) s 






Lefever et al. (2010) s 






Henrichs et al. (2009)™ au 






Catanzaro (2008)° 


51 


HD 306795 


McSwain (2008) s ™ 






Favata et al. (2009)* 






McSwain ct al. (2008)" 


19 


ALS 15211 


Gagne et al. (2011) s 


52 


HD 25558 6 


Lefever et al. (2010) sr 






Naze et al. (2012a) 6 


53 


HD 35298 " 


Landstreet et al. (2007) s 






Naze et al. (2011)* 






Bychkov et al. (2005)™ 


20 


HD 122451 " 


Ausseloos et al. (2006) s 






Yakunin et al. (2011) 6 






Alecian et al. (2011) ™ 


54 


HD 130807 


Alecian et al. (2011) s ™ 






H. Henrichs (priv. com.)" 


55 


HD 142990 


Cidale et al. (2007) a 






Raassen et al. (2005)* 






Bychkov et al. (2005)™ 


21 


HD 127381 


Henrichs et al. (2012) srbau 






Shore et al. (2004)"" 


22 


ALS 3694 


Bagnulo et al. (2006) 6 


56 


HD 37058 


Glagolevskij et al. (2007) 3 






Landstreet et al. (2007) 8 






Pedersen (1979) r 






Huang & Gies (2006) r 






Bychkov et al. (2005)™ 


23 


HD 163472 








Ramirez et al. (2004)* 






Neiner et al. (2012) r6 


57 


HD 35502 


Landstreet ct al. (2007) s 






C. Neiner (priv. com.)" 






Romanyuk & Kudryavtsev (2008) 6 






Oskinova et al. (2011)* 






Bohlender et al. (in prep)™ a 


24 


HD 96446 


Neiner et al. (2012b) sr6a 






Grillo et al. (1992)* 



Following on the next page 
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Table 2: Continued 



ID 


Star 


Ref. 


ID 


Star 


Ref. 


(1) 


(2) 


(3) 


(1) 


(2) 


(3) 






Shore & Brown (1990)" 


58 


HD 176582 


Bohlender & Monin (2011) sr,,a 


25 


HD 66765 


Cidale et al. (2007) a 

Alecian et al. (MiMeS in prep) 7,6 " 


59 


HD 189775 


Lyubimkov et al. (2002) a 
Bohlender et al. (priv. com.) ri> 


26 


HD 64740 


Bohlender & Landstreet (1990) sr 
Shore & Brown (1990)" 


60 


HD 61556 


Rivinius et al. (2003)"" 
Shultz et al. (in prep) s6a 






Peralta et al. (MiMeS in prep) ba 


61 


HD 175362 


Leone & Manfre (1997)" 






Drake et al. (1994)* 






Bychkov et al. (2005) r6 


27 


ALS 15956 


Bagnulo et al. (2006) st 
Naze et al. (2011)* 






Shore et al. (2004)°" 
Grillo et al. (1992)* 


28 


ALS 9522 


Alecian et al. (2008b) srf,[1 
Guarcello et al. (2012)* 


62 


HD 105382 


Briquet et al. (2001) sra 
Alecian et al. (2011) 6 


29 


HD 36982 " 


Wolff et al. (2004) r 


63 


HD 125823 


Bohlender et al. (2010) sr6 






Petit & Wade (2012) sf,a 


61 


HD 36526 


Landstreet et al. (2007) s 






Stelzer et al. (2005)* 






Bychkov et al. (2005) r 


30 


HD 37017" 


Bolton et al. (1998) sr 
Bohlender et al. (1987) b 
Shore & Brown (1990)" 
Oskinova et al. (2011)* 






Romanyuk & Kudryavtsev (2008) 6 


31 


HD 37479 


Hunger et al. (1989)" 
Townsend et al. (2010) r 
Oksala et al. (2012) 6a 
Shore & Brown (1990)" 
Sanz-Forcada et al. (2004)* 









s Stellar parameters, r Rotational parameters, Magnetic field parameters 
a Ha proxy, " UV proxy, * X-ray proxy 
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2.1 Sample selection 

Magnetic fields in hot stars can be detected through the 
circular polarisation induced in spectral lines by the Zee- 
man effect, using various types of instruments. The bulk of 
cooler magnetic ApBp stars were generally detected with 
first-generation photo-polarimeters, measuring for example 
the degree of polarisation in the wings of a Balmer line (e.g. 
Borra & Landstreet 1980). 

However, apart from a few strongly magnetic He-strong 
stars such as rrOriE, the bulk of hot magnetic OB stars 
were detected with second generation instruments, such as 
the low resolution (R ~ a few thousands) spectropolarime- 
try optics used in FORS 1 and 2 (VLT) and the high res- 
olution (7? ~ a few tens of thousands) spectropolarime- 
ters MUSICOS, ESPaDOnS, Narval and HARPSpol at the 
TBL, CFHT, TBL and ESO-3.6m, respectively. These two 
classes of instruments differ in that low resolution spec- 
tropolarimeters are only sensitive to the disk-integrated, 
brightness- weighted longitudinal field component, whereas 
high-resolution instruments can probe field configurations 
through the rotationally induced Doppler shifts within the 
resolved line profiles (see Donati & Landstreet 2009; Petit 
2011). 

We use the existing compilations of ApBp stars (e.g. 
Bychkov et al. 2005; Landstreet et al. 2007; Romanyuk & 
Kudryavtsev 2008) as well as an exhaustive review of the lit- 
erature to identify hot stars with confirmed field detections, 
which we complement with new detections from the MiMeS 
project. 

Some concerns have recently been raised about claimed 
magnetic detections (usually near the 3<r level) obtained 
with the FORS instruments that were not reproduced with 
other high-resolution instruments (see Silvester et al. 2009; 
Shultz et al. 2012). Bagnulo et al. (2012) performed an in- 
depth study of the complete set of FORS circular polarisa- 
tion measurements in the ESO archive, exploring the effect 
of various data reduction procedures and carefully consid- 
ering all known sources of uncertainties. Using their new 
prescription for FORS data analysis, most of the claimed 
marginal detections were found to have very low signifi- 
cance, in agreement with the results from high-resolution 
instruments. They also provided updated longitudinal field 
values and new magnetic detection statuses for stars that 
were reported magnetic in the literature at the < 6cr level. 
We therefore base our selection on these new detection sta- 
tuses for stars that were only detected with the FORS in- 
struments. 

It is worth noting that stars with chemical abundance 
peculiarities can have effective temperatures that do not re- 
flect their spectral types, as the latter is determined from 
spectral morphology. In particular He-strong/weak stars, 
which form the majority of the cooler part of our sample, 
are identified by their unusually strong/weak He lines, lines 
that are the basic means to classify B-type stars. Given that 
photometric/spectral effective temperature determinations 
are not always readily available, it is therefore difficult to 
assess the completeness of our sample at the low temper- 
ature boundary. We therefore consider all magnetic stars 
with spectral type B5 and earlier, as well as additional stars 
of later spectral type known to have effective temperatures 
greater than 16 kK. We believe the sample at these temper- 



atures (and above) to be a substantially complete represen- 
tation of the currently known hot magnetic stars. 

Although we consider a detailed review of the large sam- 
ple of stars evaluated for inclusion in Table 1 beyond the 
scope of this work, two noteworthy objects require a brief 
mention. The first of these is the Be star ui Ori, reported 
to be magnetic by Neiner et al. (2003c) based on MuSi- 
CoS longitudinal field measurements. Recently, Neiner et al. 
(2012a) have retracted this claim based on new ESPaDOnS 
and Narval measurements. The second is £ Ori A, reported 
to be magnetic by Bouret et al. (2008). While no single ob- 
servation of this star yields a significant magnetic detection, 
overall we consider the evidence presented by Bouret et al. 
(2008) to be sufficiently compelling that we retain this star 
in our list. Note that ("Ori A occupies a unique position in 
the magnetic confinement-rotation diagram (see § 3). 

2.2 Physical parameters 

Effective temperatures and surface gravities (columns 6 and 
7 of Table 1) were retrieved from the literature. An s su- 
perscript in column (6) indicates stellar parameters that 
were determined by modern spectral modelling, with NTLE 
model atmospheres such as CMFGEN, tlusty or FASTWIND 
for the hotter stars, or such as LTE ATLAS models with the 
polarised radiative transfer code zeeman for the cooler stars 
(Hillier & Miller 1998; Lanz & Hubeny 2003; Puis et al. 
2005; Kurucz 1979; Landstreet 1988; Wade et al. 2001). For 
the other stars, temperatures and gravities were generally 
derived from photometry combined with spectral type cal- 
ibrations. Some details are given in the notes of Appendix 
A in cases where significant discrepancies were found in the 
literature values or when we had to estimate log g from the 
luminosity class. 

When modern spectral modelling is available, we use the 
literature value for the luminosity, radius and mass (columns 
8, 9 and 10). The luminosity is generally obtained through 
a distance estimate and photometry, and the spectroscopic 
mass is derived from the surface gravity and radius, unless 
a better estimate is available from a binary orbit. 

For most of the remaining stars, marked with super- 
script p or c in column (8), we derive the luminosity from 
photometry (see §2.2.1) using tabulated bolometric correc- 
tions, or using the spectral energy distribution (SED) fitting 
code Chorizos (Mai'z-Apellaniz 2004) for stars with suffi- 
cient photometric data. 

In Figure 1, we locate the magnetic OB stars on the 
HR diagram. The symbol shapes represent the O-type stars 
(circles), B-type stars hotter than 22 kK (squares), those be- 
tween 22 kK and 19 kK (triangles) and those that are cooler 
than 19 kK (pentagons), and known Herbig Be stars (HeBe; 
diamonds). The luminosity classes are colour coded. The 
labels refer to the identification numbers in column (1) of 
Table 1. 

The position of the spectral types, from the calibrations 
of Martins et al. (2005) for the O-type stars and de Jager & 
Nieuwenhuijzen (1987) for the B-type stars, are indicated on 
the dark grey line that runs approximately mid-way between 
the zero-age main sequence and the terminal-age main se- 
quence; the main sequence itself is shown by the light grey 
shaded area (from the galactic evolutionary tracks of Brott 
et al. 2011). 
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Figure 1. Location of the magnetic stars in the HR diagram. The labels refer to ID sequence number listed in column (1) of Table 
1. The various symbol shapes represent effective temperature ranges and colours denote luminosity classes, as indicated in the legend. 
The shaded region shows the main sequence, from ZAMS to TAMS (from the galactic evolutionary tracks of Brott et al. 2011). The 
grey line shows the mid-way main sequence with spectral types calibrations from Martins et al. (2005) for O-type stars and de Jager & 
Nieuwenhuijzen (1987) for B-type stars. 



2.2.1 Luminosity derivation 

For each star in our sample without modern spectral 
modelling, Table 3 gathers visual magnitudes and colours 
(columns 3-5) in the Johnson UBV system, from the com- 
pilations of Mermilliod (2006) and Reed (2005) 3 . We also 
provide RJHK magnitudes (columns 6-9) from the NO- 
MAD catalog (Zacharias et al. 2005), which will be used 
below for SED fitting with Chorizos. 

For all these stars, we derive the luminosity using bolo- 
metric corrections (BC) and extinction (Av) evaluated from 
the intrinsic colour (B — V)o- The results are compiled in 
Table 4. The distance modulus (DM; column 4) is estimated 
using either Hipparcos parallax measurements or a distance 
estimate from an association with a stellar cluster. The Hip- 
parcos distances are corrected for Lutz-Kelker-type effects 

3 As no Johnson UBV measurements are available for HD 61556 
(ID 60) and HD 306795 (ID 51), we use Stromgren photometry 
from Hauck & Mermilliod (1998) with the transformation given 
by Turner (1990). 



(Lutz & Kelker 1973) using the technique described by Maiz- 
Apellaniz (2001; 2005) updated to the new reduction of the 
Hipparcos data (van Leeuwen 2007) by Mafz Apellaniz et al. 
(2008). 

The theoretical BC and (B — V)o (columns 5 and 6) 
are determined from a smooth interpolation of the grids pro- 
vided by Martins et al. (2005), and Martins & Plez (2006) 
for the O-type stars and Lanz & Hubeny (2007) for the 
B-type stars. We use an extinction i?v = 3.1 to derive 
the extinction Av = Rv E(B — V) (column 7). The ab- 
solute visual magnitude (My = V — Av — DM), the bolo- 
metric magnitude (Mboi = My + BC) and the luminosity 
[log(L*/L Q ) = (Mboi,© — Mboi)/2.5] are given in columns 8 
to 10. 

With a typical uncertainty of 2 000 K in T e g and 0.3 dex 
in log <j, we estimate an uncertainty of 0.2 and 0.02 mag in 
BC and (B — V)o, respectively. Given the wide range of Rv 
often encountered in the literature for OB stars, we adopt 
a conservative error in Av of 0.25 mag. In most cases, BC, 
Ay and DM contribute equally to the uncertainty, leading 
to 0.2-0.3 dex for the luminosity. In five cases (ID: 36, 37, 40, 
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Table 3. Photometry of magnetic stars without modern spectral modelling (§ 2.2.1). 



ID 


Star 


V 


(B-V) 


(17 — B) 


R 


J 


H 


K 






±0.01 mag 


±0.01 mag 


mag 


±0.02 mag 


±0.02 mag 


±0.02 mag 


±0.02 mag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


20 


HD 122451 1 


0.61 


-0.24 












22 


ALS 3694 


10.38 


0.33 


-0.55 ±0.02 


10.26 


9.62 


9.49 


9.47 


25 


HD 66765 


6.62 


-0.16 


-0.8 ±0.2 


6.69 


6.89 


7.01 


7.01 


27 


ALS 15956 


10.80 


0.05 


-0.69 ±0.02 


10.70 


10.62 


10.64 


10.66 


29 


HD 36982 


8.45 


0.13 


-0.57±0.02 


8.42 


7.74 


7.64 


7.47 


31 


HD 37479 


6.66 


-0.18 


-0.87±0.01 


6.40 


6.97 


6.95 


6.95 


32 


HD 149277 1 


8.38 


0.04 


-0.658 ±0.009 


8.43 


8.26 


8.26 


8.28 


35 


HD 136504 1 


3.37 


-0.18 


-0.74±0.02 


3.44 


3.99* 


3.93* 


4.13* 


36 


HD 156424 


8.72 


0.06 












37 


HD 156324 


8.75 


0.10 












38 


HD 121743 


3.82 


-0.22 












40 


HD 186205 t 


8.53 


0.05 


-0.60±0.02 


8.50 


8.40 


8.48 


8.48 


41 


HD 67621 


6.33 


-0.20 


-0.802 ±0.004 


6.41 


6.73 


6.86 


6.84 


43 


HD 35912 


6.41 


-0.18 












44 


HD 66522 


7.20 


0.05 


-0.624 ±0.005 


7.19 


7.03 


7.02 


7.01 


48 


HD 58260 t 


6.74 


-0.13 


-0.756 ±0.006 


6.79 


6.98 


7.11 


7.12 


49 


HD 36485 t 


6.85 


-0.16 


-0.72 ±0.01 


6.92 


7.17 


7.23 


7.28 


50 


HD 208057 1 


5.07 


-0.17 


-0.7±0.1 


5.16 


5.39 


5.49 


5.54 


51 


HD 306795 


10.62 


0.02 












52 


HD 25558 t 


5.32 


-0.08 


-0.57±0.01 


5.35 


5.48 


5.58 


5.53 


53 


HD 35298 + 


7.89 


-0.14 


-0.6 ±0.2 


7.94 


8.12 


8.17 


8.27 


54 


HD 130807 


4.32 


-0.15 


-0.62 ±0.01 


4.38 


4.78* 


4.72 


4.73 


55 


HD 142990 1 


5.43 


-0.09 


-0.653 ±0.003 


5.46 


5.58 


5.67 


5.65 


56 


HD 37058 t 


7.30 


-0.13 


-0.79 ±0.02 


7.23 


7.60 


7.73 


7.75 


57 


HD 35502 t 


7.35 


-0.04 


-0.54±0.04 


7.32 


7.39 


7.42 


7.43 


59 


HD 189775 


6.14 


-0.19 


-0.66 ±0.02 


6.20 


6.46 


6.53 


6.64 


60 


HD 61556 t 


4.78 


-0.14 












61 


HD 175362 1 


5.37 


-0.15 


-0.7±0.1 


5.43 


5.62 


5.66 


5.68 


62 


HD 105382 t 


4.46 


-0.16 


-0.68 ±0.01 


4.51 


5.12 


4.95 


4.87 


63 


HD 125823 


4.37 


-0.19 












64 


HD 36526 


8.31 


-0.11 













* Because of their brightness, these stars have uncertainties of ±0.24 mag according to the 2mass specifications. 



53 and 57) the luminosity error estimate from the bolomet- 
ric correction method is more than 0.4 dex, given the large 
uncertainty in distance. 

For the stars with a complete set of UBVRJHK pho- 
tometry, we perform SED fitting using the Bayesian (spec- 
tro)photometric code Chorizos. The results are presented 
in Table 5. In the latest Chorizos version, the user can 
select distance to be an independent parameter, by apply- 
ing atmosphere models (tlusty for OB stars) calibrated in 
luminosity with the help of Geneva stellar evolution tracks 
(excluding rotational effects). The parameters of such mod- 
els are the logarithmic distance (column 3), the extinction 
(here fixed at -R5495 = 3.1 4 ), the reddening (transformed to 
Ay in column 4), the effective temperature (here fixed to the 
literature estimate) and the luminosity class. The distance 
prior probability range has been left relatively wide around 
the Hipparcos or cluster-estimated value and the luminos- 
ity class prior probability was based on the gravity estimates 
used for the bolometric correction approach (column 3 of ta- 

4 The extinction law is defined by the monochromatic quantity 
-R5495 = A54g5/_B(4405 — 5495) instead of a band-integrated one 
such as Ry = Ay I E(B — V), because the former depends only 
on the properties of the dust while the latter also depends on the 
input SED and the amount of dust present along the line of sight. 
See Mafz Apellaniz (2012) for details. 



ble 4) with an interval of 1.5a. From these fitted values, we 
can derive an estimate of the surface gravity (column 5), the 
evolutionary mass (column 6) and the luminosity (column 
7). 

Good fits to the photometry are achieved for the 13 
stars displayed in Table 5, leading to better estimates of their 
luminosity (especially for the two stars with the largest un- 
certainty with the bolometric correction approach). Poorer 
fits were obtained for the remaining 9 stars with complete 
photometry. Incompatibility between optical and near-IR 
photometry, probably due to near-IR excess, could be a pos- 
sible cause of the discrepancy. Therefore, for the stars with 
good fits, we use the luminosity, gravity and mass derived 
from Chorizos. For the remaining stars, we opt for the bolo- 
metric correction luminosity determination. 

2.3 Rotational and magnetic parameters 

Monitoring of the disk-integrated longitudinal field varia- 
tions provides a natural and direct way to determine ro- 
tational periods for magnetic stars (in the context of the 
Oblique Rotator Model; e.g. Stibbs 1950). Photometric and 
spectral variability associated with the magnetic field also 
provide a convenient and easy way to determine periods, 
even though some ambiguity can exist between e.g. short 
rotational periods and long pulsation periods. 
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Table 4. Luminosity determination based on bolometric correction and extinction from intrinsic colours (§ 2.2.1). 



ID Star logg DM BC (B — V)o A v M v M bol log(L*/L ) 

cgs mag ±0.2 mag ±0.02 mag ±0.25 mag mag mag 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



20 


HD 122451 1 


3.5 ±0.4 




-2.48 






-3.8±0.5 a 


-6.3 ±0.5 


4.4 ±0.2 


22 


ALS 3694 


4.0±0.4 


11.2±0.6 C 


-1.98 


-0.18 


1.57 


-2.4 ±0.7 


-4.4 ±0.8 


3.7±0.3 


25 


HD 66765 


4.0±0.2 


8.4 ±0.4 


-1.98 


-0.18 


0.06 


-1.9±0.5 


-3.8±0.6 


3.4±0.2 


27 


ALS 15956 


4.0±0.4 


12.8±0.6 C 


-2.32 


-0.20 


0.78 


-2.8 ±0.7 


-5.1±0.7 


3.9 ±0.3 


29 


HD 36982 


4.0±0.2 


8.3±0.2 C 


-2.21 


-0.19 


1.01 


-0.8±0.3 


-3.0±0.4 


3.1 ±0.2 


31 


HD 37479 


4.0±0.2 


8.5 ± 0.4 c 


-2.31 


-0.20 


0.07 


-1.9±0.5 


-4.2 ±0.5 


3.6 ±0.2 


.32 


HD 149277 + 


4.0±0.4 


10.7±0.6 C 


-2.21 


-0.19 


0.72 


-3.1 ±0.7 


-5.3±0.7 


4.0 ±0.3 


35 


HD 136504 1 


4.0±0.2 


6.1 ±0.3 


-2.21 


-0.19 


0.04 


-2.7 ±0.4 


-5.0±0.4 


3.9 ±0.2 


.36 


HD 156424 


4.0±0.3 


10±l c 


-2.21 


-0.19 


0.78 


-2±1 


-4±1 


3.7±0.4 


.37 


HD 156324 


4.0±0.3 


10±l c 


-2.21 


-0.19 


0.92 


-2±1 


-5±1 


3.7±0.4 


38 


HD 121743 


4.0±0.3 


6.03 ±0.07 


-2.21 


-0.19 


0.00 


-2.2 ±0.3 


-4.4 ±0.4 


3.7±0.2 


10 


HD 186205 t 


4.0±0.2 


14 ±1 


-1.98 


-0.18 


0.70 


-6±1 


-8±1 


5.1 ±0.5 


41 


HD 67621 


4.0±0.3 


8.0±0.2 


-1.85 


-0.17 


0.00 


-1.7±0.3 


-3.5±0.5 


3.3±0.2 


43 


HD 35912 


4.0±0.1 


8.2 ±0.7 


-1.78 


-0.16 


0.00 


-1.8±0.7 


-3.5±0.7 


3.3±0.3 


M 


HD 66522 


4.0±0.4 


8.3 ±0.4 


-1.71 


-0.16 


0.64 


-1.8 ±0.4 


-3.5±0.5 


3.3 ±0.2 


48 


HD 58260 t 


3.8±0.3 


9.7±0.5 


-1.97 


-0.18 


0.15 


-3.2 ±0.6 


-5.1 ±0.6 


3.9 ±0.3 


49 


HD 36485 t 


4.2±0.2 


8.4±0.7 C 


-1.99 


-0.18 


0.06 


-1.6±0.7 


-3.6±0.8 


3.3±0.3 


50 


HD 208057 + 


3.9±0.2 


6.4±0.1 


-1.85 


-0.17 


0.00 


-1.4±0.3 


-3.2±0.5 


3.2±0.2 


51 


HD 306795 


3.9±0.2 


11.6±0.6 C 


-1.71 


-0.16 


0.55 


-1.6 ±0.7 


-3.2±0.7 


3.2±0.3 


52 


HD 25558 + 


4.0±0.2 


6.5 ±0.1 


-1.64 


-0.15 


0.21 


-1.4±0.3 


-3.0±0.4 


3.1 ±0.2 


5.3 


HD 35298 + 


4.0±0.3 


12 ±1 


-1.41 


-0.14 


0.00 


-4±1 


-5±1 


3.9 ±0.5 


54 


HD 130807 


4.2±0.2 


5.5 ±0.2 


-1.72 


-0.15 


0.00 


-1.2±0.3 


-2.9 ±0.4 


3.1 ±0.2 


55 


HD 142990 t 


4.2±0.2 


6.16±0.09 


-1.64 


-0.15 


0.17 


-0.9±0.3 


-2.6±0.4 


2.9 ±0.2 


56 


HD 37058 + 


3.8±0.2 


8.5 ± 0.4 c 


-1.71 


-0.16 


0.09 


-1.3±0.5 


-3.0±0.6 


3.1 ±0.2 


57 


HD 35502 t 


4.0±0.3 


9±2 


-1.41 


-0.14 


0.30 


-2 ±2 


-4 ±2 


3.3±0.8 


59 


HD 189775 


4.0±0.3 


6.9±0.1 


-1.41 


-0.14 


0.00 


-0.8±0.3 


-2.2±0.4 


2.8±0.1 


60 


HD 61556 t 


4.0±0.3 


5.2 ±0.2 


-1.24 


-0.12 


0.00 


-0.4±0.3 


-1.7±0.4 


2.6±0.1 


61 


HD 175362 t 


4.0±0.3 


5.60 ±0.08 


-1.24 


-0.12 


0.00 


-0.2 ±0.3 


-1.5 ±0.4 


2.5 ±0.1 


62 


HD 105382 t 


4.0±0.2 


5.7±0.2 


-1.64 


-0.15 


0.00 


-1.2±0.3 


-2.9 ±0.4 


3.0 ±0.2 


6.3 


HD 125823 


4.0±0.2 


5.73 ±0.06 


-1.85 


-0.17 


0.00 


-1.4±0.3 


-3.2±0.4 


3.2±0.1 


64 


HD 36526 


4.0±0.3 


8.4±0.7 C 


-1.41 


-0.14 


0.08 


-0.1 ±0.7 


-1.5±0.8 


2.5 ±0.3 



a Prom the SB2 analysis of Ausseloos et al. (2006). 

c Distance estimates from associations with stellar clusters (Hipparcos otherwise). 



Table 5. Luminosity determination based on SED fitting with CHORIZOS (§ 2.2.1). 



ID 


Star 


DM 


A v 


logg 


M* 


log(L*/L ) 






mag 


mag 


cgs 


M e 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


25 


HD 66765 


8.6 ±0.5 


0.21 ±0.02 


3.9 ±0.2 


7.5 ±0.6 


3.6 ±0.2 


27 


ALS 15956 


13.6 ±0.6 


0.92 ±0.02 


3.6 ±0.2 


U±l 


4.3 ±0.2 


35 


HD 136504 1 


5.6 ±0.5 


0.19 ±0.04 


4.0 ±0.2 


8.6±0.7 


3.8 ±0.2 


40 


HD 186205 t 


9.9 ±0.5 


0.75 ±0.02 


4.0 ±0.2 


7.4±0.6 


3.5 ±0.2 


48 


HD 58260 1 


10.0 ±0.5 


0.24 ±0.01 


3.5±0.2 


9±1 


4.1 ±0.2 


49 


HD 36485 f 


8.6 ±0.4 


0.19 ±0.02 


4.0±0.1 


7.1 ±0.4 


3.5 ±0.1 


50 


HD 208057 f 


7.2 ±0.5 


0.12 ±0.02 


3.9 ±0.2 


7.1 ±0.6 


3.6 ±0.2 


53 


HD 35298 f 


9.7±0.6 


0.12 ±0.02 


3.8 ±0.2 


5.6±0.5 


3.2 ±0.2 


54 


HD 130807 


5.3 ±0.3 


0.13±0.01 


4.1 ±0.1 


5.7±0.2 


3.1±0.1 


56 


HD 37058 1 


9.4±0.5 


0.07 ±0.03 


3.8 ±0.2 


6.6±0.5 


3.5 ±0.2 


57 


HD 35502 1 


8.9 ±0.6 


0.40 ±0.01 


3.8 ±0.2 


5.7±0.5 


3.3 ±0.2 


59 


HD 189775 


7.9 ±0.5 


0.00 ±0.02 


3.8 ±0.2 


5.5±0.5 


3.2 ±0.2 


61 


HD 175362 t 


7.2 ±0.4 


0.08 ±0.01 


3.7±0.2 


5.3 ±0.4 


3.2±0.1 



In Table 1, column (11) gives the rotational period in 
days. When no period is available we use the measured v sin i 
(column 12) as a lower limit to the equatorial velocity. In 
two cases (ALS 15211, ID 19; ALS 15956, ID 27), no vsini 
measurements are available, due to a lack of high-resolution 



spectra. These stars would be prime candidates for further 
monitoring. 

In four cases (HD 96446, HD 136504, HD 58260 and 
HD 37058; ID 24, 35, 48 and 56) more than one period is 
reported in the literature. In these cases, we use the longest 
period for a lower limit on the equatorial velocity. We pro- 
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vide the magneto-spheric calculations for the alternative pe- 
riods in the notes of Appendix A. 

Column (13) gives the estimated polar strength (B p ) of 
the surface dipole in kilogauss. When only longitudinal mag- 
netic measurements are available, we use a value of three 
times the strongest longitudinal field measurement (corre- 
sponding to a conservative limb-darkening coefficient of 
0.6), setting a lower limit on the dipolar field strength. A 
superscript m in column (13) indicates stars that are known 
to possess a magnetic field with a significant contribution 
from multipole components higher than a simple dipole. 



3 TWO-PARAMETER CLASSIFICATION OF 
MAGNETOSPHERES 

3.1 Alfven radius Ra vs. Kepler co-rotation 
radius Rk 

The high luminosity of massive stars drives powerful, high- 
speed stellar winds. MHD simulation studies (e.g. ud-Doula 
& Owocki 2002; ud-Doula et al. 2008) show that the overall 
net effect of a large-scale, dipole magnetic field in divert- 
ing such a wind can be well characterised by a single wind 
magnetic confinement parameter, 



B eq R, 
Mb=o 



(1) 



where -B eq = B p /2 is the field strength at the magnetic equa- 
torial surface radius _R„, and Mb=o and Voo are the fiducial 
mass-loss rate and terminal speed that the star would have 
in the absence of any magnetic field. 

This confinement parameter sets the scaling for the ra- 
tio of the magnetic to wind kinetic energy density. For a 
dipole field, the r -6 radial decline of magnetic energy den- 
sity is much steeper than the r~ 2 decline of the wind's mass 
and energy density; this means the wind always dominates 
beyond the Alfven radius Ra (ud-Doula et al. 2008), given 
by the approximate general scaling, 



-Ra 
R, 



0.3 + f> + 0.25) 



1/4 



(2) 



Magnetic loops extending above Ra are drawn open by the 
wind, while those with an apex below Ra remain closed. 
Indeed, the trapping of wind upfiow from opposite foot- 
points of closed magnetic loops leads to strong collisions that 
may form X-ray emitting, magnetically confined wind shocks 
(MCWS Babel & Montmerle 1997a,b, see §5.2). In models 
with negligible rotation, the post-shock material eventually 
cools and falls back onto the star, leading to a relatively 
complex, dynamic pattern of infall and wind outflow (see 
e.g. lower row of Figure 9 of ud-Doula et al. 2008, also Fig- 
ure 2). 

For the simple 2D axisymmetric case of a magnetic 
dipole that is aligned with a star's rotation axis, ud-Doula 
et al. (2008) extended these MHD simulation studies to ex- 
plore the additional effect of stellar rotation. They found it 



5 In the case of a dipolar field, the dipole strength can be ex- 
pressed as B p > 4^=^|<B Z )| max, where e is the limb-darkening 
coefficient and |(-B z )| max is the maximum of the disk-integrated 
longitudinal field variation (Preston 1967). 



convenient to cast results in terms of the ratio of the rota- 
tion speed V I0 t to orbital speed V or b at the equatorial surface 
radius 



W 



Vrot 
Voib 



U-R, 



s/GM*/R, 



(3) 



where the latter equality expresses this ratio in terms of 
the angular rotation frequency u), with M* the stellar mass. 
To avoid the complications associated with a rotationally 
distorted, oblate stellar surface, ud-Doula et al. (2008) re- 
stricted their simulations to cases with W ^ 0.5. But if we 
associate -R* with the actual equatorial radius for the given 
rotation rate uj , then even for more rapid, near-critical rota- 
tion, W simply compares the star's equatorial rotation speed 
to the speed V Q ib needed to reach Keplerian orbit near this 
equatorial surface 6 . 

In a magnetic star, torques from the magnetic field 
on any wind outflow can maintain rigid-body co-rotation 
up to roughly the Alfven radius, so that the azimuthal 
speed of the confined wind plasma increases with radius 
as Vj, = Vrotr/R*. The outward centrifugal force from such 
rigid-body rotation will balance the inward force of gravity 
at the Kepler corotation radius, 



Rk 



1/3 



w 



-2/3 



(4) 



Together the two parameters 77* and W thus define the rel- 
ative locations of the Alfven and Kepler radii with respect 
to the equatorial radius. 



3.2 Dynamical vs. Centrifugal Magnetospheres 

For the simple case of field-aligned rotation, ud-Doula et al. 
(2008) carried out an extensive MHD simulation parameter 
study varying both W and 77*. For n, < 1, the field exerts 
only a modest perturbation on the wind; but for n* > 1, out- 
flow near the magnetic equator is trapped within the Alfven 
radius by closed magnetic loops, forming a wind-fed circum- 
stellar magnetosphere. It was found that the dynamical evo- 
lution of this trapped magnetospheric material depends cru- 
cially on the rotation parameter W, and specifically on the 
relative magnitude of the associated Kepler vs. the Alfven 
radii. 

In a simplified, schematic form, Figure 2 here illustrates 
that, depending on the relative positions of Rk vs. Ra, re- 
gions of trapped equatorial material can be alternatively 
characterised as forming a dynamical vs. centrifugal mag- 
netosphere (DM vs. CM). As sketched in the upper panel of 
Figure 2, for slowly rotating stars with Ra < Rk, material 
trapped in closed magnetic loops falls back to the star on a 
dynamical timescale, forming a DM (Sundqvist et al. 2012). 
In contrast, the lower panel shows that, for the more rapidly 
rotating case with Ra > -Rk, material caught in the region 



6 For critical rotation (W = 1), R t = 3R p /2, where R p is the 
fixed polar radius. In terms of the associated critical rotational fre- 
quency u; cr it = ^J&GM/ 27 R3 , one can alternatively define a crit- 
ical rotation ratio f! = u]/ui cr - lt . We then find W = Q(2R*/3R P ), 
with R*/R p = 3cos[(cos _1 [fi] +7r)/3]/fi (Collins & Harrington 
1966). 
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Figure 2. Sketch of the regimes for dynamical vs. centrifugal 
magnetospheres (DM vs. CM). The top panel illustrates the case 
of a slowly rotating star with Kepler radius beyond the Alfvcn 
radius (Rk > Ra)\ the lack of centrifugal support means that 
trapped material falls back to the star on a dynamical timescale, 
forming a DM, with colour illustrating the rough time-averaged 
distribution of density. The lower panel is for a star with more 
rapid rotation and Rk < Ra, leading then to a region between 
these radii where a net outward centrifugal force against gravity 
is balanced by the magnetic tension of closed loops; this allows 
material to build up to the much higher density of CM. 

between Ra and Rk is centrifugally supported against in- 
fall, and so builds up to a much denser CM (for a given 
fiducial mass- loss rate). Even for such rapid rotators, the 
inner regions below Rk again have the infall of a DM, but 
the plasma density, and thus any circumstellar emission, is 
much lower than that of the CM region. 

Indeed, since the much longer confinement time allows 
material to accumulate to high density even if the feeding 
by the wind mass flux is weak, such a CM can exhibit rota- 
tionally modulated line emission even in the relatively low- 
luminosity, but strongly magnetic Bp stars, so long as the 
stellar rotation is sufficient to give Rk < Ra (Townsend 
& Owocki 2005; Townsend et al. 2005). For slowly rotating 
magnetic stars with a DM, accumulating sufficiently high 
density plasma for line emission requires a much stronger 
wind to overcome the dynamical timescale leakage of infall 



back onto the star. For the luminous, slowly rotating mag- 
netic O-type star HD 191612 (ID 4), Sundqvist et al. (2012) 
showed that the emission from its wind-fed DM matches its 
observed Ha emission quite well. 

The transition from stars with a pure DM to those with 
a CM occurs near Rk = Ra; from equations (2) and (4), the 
associated transition value W t for the rotation fraction is 

W t = [o.3 + (tj. + 0.25) 1/4 ] ~ A ' 2 , (5) 
which in the strong confinement limit, r]» 1, simply re- 

Tir -3/8 

quires Wt ~ V* 

Figure 3 plots our sample of magnetic stars in the mag- 
netic confinement-rotation diagram, a log-log plane with 
Rk/ R* increasing downward on the ordinate vs. Ra/ R* in- 
creasing to the right along the abscissa. As detailed in the 
next subsection (§ 3.3), the placement of the individual stars 
depends on inference of the relevant parameters that set the 
magnetic confinement r/* (noted on the top axis) and ro- 
tation fraction W (on the right axis). The vertical line at 
77* = 1 (Ra/R* ~ 1.3) separates weakly magnetised winds 
at the far left from the broad domain of stars with significant 
magnetospheres, with the diagonal line separating the stars 
with a CM to the upper right from those with just a DM to 
the lower left. As detailed in § 3.4, the additional upper and 
right axes refer to associated stellar spindown properties, 
namely the stellar spindown timescale (rj) and the maxi- 
mum spindown age (t s ,max), respectively. Stars above the 
horizontal dotted line have a maximum spindown age i Simax 
that is less than one spindown time tj. 

3.3 Calculation of magnetospheric parameters 

In this section, we determine the magnetospheric parame- 
ters described in §3.1 for all the stars in the sample. Table 
6 compiles our calculations of 77*, Ra/R*, W and Rk/R* 
(equations 1 to 4) in columns (4) to (7). 

3.3.1 Wind momentum 

To compute 77* and Ra from equations 1 and 2, we need, in 
addition to the stellar radius and surface magnetic field, es- 
timates of the wind mass-loss rate Mb=o and terminal speed 
Voo- Simulation models define the confinement in terms of 
wind properties a star would have if it had no magnetic 
field. Therefore instead of making empirical estimates of the 
wind properties of each magnetic star (which are in any case 
difficult to obtain, see Sundqvist et al. 2012; Grunhut et al. 
2012c), we derive theoretical values based on inferred stellar 
parameters applied to radiation line-driven wind theory. 

Following standard theory, we take the wind terminal 
speed Voo to scale with the star's effective surface escape 
speed, 

K sc ,( 2GM f^) 1/ \ (6) 

where r e = K c L/47rGAf»c is the Eddington parameter for 
electron scattering opacity K e . For the order-unity ratio 
Kjc/Vcsc we use the factors recommended by Vink et al. 
(2000, 2001) (based on the empirical study of Lamers et al. 
1995), which declines abruptly from 2.6 to 1.3 from the 
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Figure 3. Location of magnetic massive stars in the magnetic confinement-rotation diagram; a log-log plot with Kepler radius i?K 
increasing downward and Alfven radius i?A increasing to the right. The right and upper axes give respectively the corresponding rotation 
fraction W and magnetic confinement parameter 77*. The solid lines separates the magnctosphcrc domains of weakly magnetised winds 
(with »y* < 1), dynamical magnetospheres with < Rk (DM), and centrifugal magnetospheres with i?A > (CM), as defined 
in §3.2. The additional upper and right axes give respectively the corresponding spindown timescale Tj (normalised by the value in a 
non-magnetised wind) and maximum spindown age t s ,max (normalised by the spindown time and therefore the number of spindown 
e-folds), as defined in §3.4. Stars above the dashed line have a maximum spindown age less than one spindown time. As in Figure 1, 
the symbol shapes denote spectral type, and the numbers correspond to the ID in column (1) of Table 1. The three downward arrows 
indicate two stars (ALS 15211, ID 19; ALS 15956, ID 27) for which no vsini measurement is available (e.g. W > 0), and HD 108 (ID 7) 
for which R K ~ 500H* . 



Tabic 6: List of magnetospheric parameters, magnetic spindown properties, 
and Ha, UV and X-ray proxies. 



ID 


Star 


Remark 


W 


Rk/R* 


V* 


Ra/R* 




^s,max 


\og(R A /R K ) 


Ha 


UV 


log(L x /L*) 
















Myr 


Myr 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


HD 148937 




1.1 c-1 


4.3 


4.9 


1.8 


0.82 


1.8 


-0.37 


var em 


var 


-6.6 


2 


CPD-28 2561 




1.2 e-2 


19 


> 6.4el 


> 3.1 


< 0.92 


< 4.1 


> -0.79 


var em 






3 


HD 37022 t 


SB1 


3.5 e-2 


9.4 


1.8cl 


2.4 


3.1 


10 


-0.60 


var em 


per 


-5.6 


4 


HD 191612 t 


SB2 


2.3 e-3 


57 


1.3 c2 


3.7 


0.38 


2.3 


-1.20 


var em 


var 


-6.0 


5 


NGC 1624-2 




3.8 c-3 


41 


> 1.5c4 


> 11 


< 0.24 


< 1.4 


> -0.56 


var em 




-6.6 


6 


HD 47129 t 


SB2 


> 3.1 c-1 


< 2.2 


> 6.8 c2 


> 5.4 


< 4.1 


< 4.8 


> 0.39 


em 


var 


-5.6 


7 


HD108 




8.3 e-5 


526 


> 3.5 


> 1.7 


< 0.90 


< 8.5 


> -2.50 


var em 


var 


-6.2 


8 


ALS 15218 t 




> 3.2 e-2 


< 9.9 


> 1.2c2 


> 3.6 


< 3.1 


< 10 


> -0.44 


em 




-6.3 


9 


HD 57682 




8.1 e-3 


24 


1.3 e2 


3.7 


2.2 


10 


-0.83 


var em 


var 


-6.3 


10 


HD 37742 t 


SB2 


3.3 e-1 


2.1 


2.0e-l 


1.1 


3.9 


4.4 


-0.28 


var em 


abs 


-6.7 


11 


HD 149438 




1.1 e-2 


20 


4.6 


1.8 


21 


97 


-1.05 


var abs 


per 


-6.3 


12 


HD 37061 1 


SB2 


> 2.8 e-1 


< 2.3 


1.1 c2 


3.5 


27 


< 35 


> 0.18 


abs 


stab abs 


-6.9 


13 


HD 63425 




> 1.5 e-2 


< 16 


6.3cl 


3.1 


22 


< 93 


> -0.73 


stab abs 


var 




14 


HD 66665 




2.4 e-2 


12 


1.9c2 


4.0 


12 


46 


-0.48 


stab abs 


var 




15 


HD 46328 


/3Cep 


2.3 e-1 


2.7 


> 6.2 c2 


> 5.3 


< 1.6 


< 2.3 


> 0.30 


em 


em 


-6.8 


16 


ALS 8988 


HeBe 


> 3.3 e-2 


< 9.6 


> 2.7c3 


> 7.5 


< 23 


< 78 


> -0.11 


HeBe 




-8.8 


17 


HD 47777 


HeBe 


> 1.1 c-1 


< 4.3 


> 4.8 c3 


> 8.6 


< 8.9 


< 19 


> 0.30 


HeBe 




-6.9 


18 


HD 205021 1 


SB2,/3Cep 


5.1 e-2 


7.3 


1.8 e2 


4.0 


43 


128 


-0.26 


abs 


per 


-7.2 


19 


ALS 15211 1 








> 3.0e3 


> 7.7 


< 17 




> -2.29 






-7.8 


20 


HD 122451 1 


SB2,/3Ccp 


> 1.7 e-1 


< 3.2 


> 6.5 el 


> 3.1 


< 14 


< 26 


> -0.02 


abs 


stab abs 


-7.3 


21 


HD 127381 




1.4 c-1 


3.8 


2.7c3 


7.5 


127 


254 


0.30 


stab abs 


var 




22 


ALS 3694 




> 1.3 c-1 


< 3.8 


> 1.2c5 


> 18 


< 3.0 


< 6.1 


> 0.69 








23 


HD 163472 


/3Cep 


8.4 e-2 


5.2 


5.9 e2 


5.2 


157 


390 


0.00 


abs 


per 


-8.6 


24 


HD 96446 1 




6.9 e-2 


6.0 


7.9 c4 


17 


2.2 


5.8 


0.46 


stab abs 


em 




25 


HD 66765 




3.1 e-1 


2.2 


> 1.2c4 


> 10 


< 4.8 


< 5.6 


> 0.69 


stab abs 






26 


HD 64740 




3.8 e-1 


1.9 


8.2 c5 


30 


1.7 


1.6 


1.20 


var em 


per 


-7.3 


27 


ALS 15956 








> 8.8 c3 


> 9 


< 7.4 




> -2.14 






-6.8 


28 


ALS 9522 


HeBe 


> 3.5 e-1 


< 2.0 


> 1.3 c4 


> 11 


< 1.3 


< 1.4 


> 0.74 


HeBe 




-7.3 


29 


HD 36982 




> 1.9c-l 


< 3.0 


1.3 c3 


6.3 


8.6 


< 14 


> 0.32 


stab abs 


stab abs 


-7.8 


30 


HD 37017 1 


SB2 


3.7 e-1 


1.9 


> 1.1 c5 


> 18 


< 3.7 


< 3.6 


> 0.98 


var em 


per 


-7.8 


31 


HD 37479 




3.4 e-1 


2.1 


8.9 c5 


31 


4.6 


5.0 


1.18 


var em 


per 


-5.9 


32 


HD 149277 1 


SB2 


> 2.2 e-2 


< 12 


> 2.7 e4 


> 13 


< 2.6 


< 10 


> 0.01 


abs 






33 


HD 184927 




4.7 e-2 


7.7 


1.7c4 


11 


1.7 


5.2 


0.18 


stab abs 


per 




34 


HD 37776 t 




2.4 e-1 


2.6 


3.1 e5 


23 


0.68 


0.97 


0.96 


var em 


per 




35 


HD 136504 1 


SB2,/3Cep 


> 7.3 e-2 


< 5.7 


> 4.2 c2 


> 4.8 


< 12 


< 32 


> -0.08 


abs 


abs 




36 


HD 156424 




> 2.6 e-2 


< 11 


> 5.7e2 


> 5.2 


< 15 


< 55 


> -0.34 


var em 






37 


HD 156324 




> l.Oc-1 


< 4.6 


> 4.2 e3 


> 8.4 


< 5.9 


< 13 


> 0.26 


var em 






38 


HD 121743 


/3Cep 


> 1.4 c-1 


< 3.7 


> 3.8 c2 


> 4.7 


< 18 


< 35 


> 0.10 


abs 


abs 


-7.2 


39 


HD 3360 


SPB 


1.1 c-1 


4.4 


> 2.1 c2 


> 4.1 


< 19 


< 43 


> -0.03 


abs 


per 


-7.8 


40 


HD 186205 t 




> 9.1 e-3 


< 22 


> 8.5 e3 


> 9.9 


< 6.9 


< 32 


> -0.36 


stab abs 






41 


HD 67621 




1.1 c-1 


4.4 


> 4.9 c3 


> 8.6 


< 22 


< 49 


> 0.29 


stab abs 


abs 




42 


HD 200775 f 


SB2, HeBe 


2.8 e-1 


2.3 


3.3 c3 


7.9 


3.6 


4.6 


0.53 


HeBe 




-6.3 


43 


HD 35912 




4.5 e-1 


1.7 


> 2.8 e5 


> 23 


< 4.1 


< 3.3 


> 1.13 


abs 
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Table 6: Continued 



ID 


Star 


Remark 


W 


Rk/R* 


V* 


Ra/R* 


TJ 


is, max 


\og(R A /R K ) 


Ha 


UV 


log(L x /-L*) 
















Myr 


Myr 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


44 


HD 66522 




> 1.8 e-2 


< 14 


8.1 e3 


9.8 


30 


< 124 


> -0.18 


abs 






45 


HD 182180 




5.8 e-1 


1.4 


2.9 c6 


41 


4.5 


2.5 


1.46 


var em 


abs 




46 


HD 55522 




1.1 c-1 


4.4 


> 1.4c5 


> 19 


< 21 


< 48 


> 0.64 


abs 






47 


HD 142184 




5.0 e-1 


1.6 


4.3 c6 


45 


9.5 


6.7 


1.46 


var em 




-6.7 


48 


HD 58260 1 




> 2.7 e-2 


< 11 


> 5.2 e4 


> 15 


< 0.29 


< 1.1 


> 0.14 


abs 


cm 




49 


HD 36485 1 


SB2 


2.7 e-1 


2.4 


3.5 e5 


24 


1.6 


2.1 


1.01 


var em 


var 




50 


HD 208057 1 


SPB 


3.8 e-1 


1.9 


> 9.6 c2 


> 5.9 


< 18 


< 18 


> 0.49 


abs 


abs 




51 


HD 306795 




> 1.5 c-1 


< 3.6 


> 1.9c5 


> 21 


< 2.8 


< 5.5 


> 0.77 


abs 






52 


HD 25558 t 


SPB 


> 5.4 e-2 


< 7.0 


> 3.1 e2 


> 4.5 


< 185 


< 542 


> -0.19 








53 


HD 35298 + 




3.2 c-1 


2.1 


> 1.4c6 


> 34 


< 1.7 


< 1.9 


> 1.22 


stab abs 






54 


HD 130807 




> 4.4 e-2 


< 8.0 


> 5.3 e4 


> 15 


< 20 


< 64 


> 0.29 








55 


HD 142990 1 




2.8 e-1 


2.4 


> 1.2c6 


> 33 


< 9.3 


< 12 


> 1.15 


var em 


per 




56 


HD 37058 1 




3.9 e-2 


8.6 


> 6.2 c4 


> 16 


< 3.7 


< 11 


> 0.27 








57 


HD 35502 1 


SB2 


7.4 e-1 


1.2 


> 7.5 e5 


> 29 


< 2.1 


< 0.61 


> 1.39 


var em 




-6.0 


58 


HD 176582 




2.3 e-1 


2.7 


1.7c6 


36 


6.7 


9.9 


1.14 


var em 






59 


HD 189775 




2.2 e-1 


2.7 


> 3.8 c5 


> 25 


< 3.8 


< 5.7 


> 0.96 








60 


HD 61556 t 




1.7 e-1 


3.3 


1.1 c6 


32 


14 


26 


1.00 


stab abs 






61 


HD 175362 t 




1.9 e-1 


3.0 


> 1.2c7 


> 59 


< 0.74 


< 1.2 


> 1.29 


stab abs 


var 


-6.3 


62 


HD 105382 t 




2.8 e-1 


2.3 


7.4 c4 


16 


13 


16 


0.86 


abs 




-6.8 


63 


HD 125823 




4.1 e-2 


8.4 


> 1.0e4 


> 10 


< 14 


< 47 


> 0.09 


var abs 


per 




64 


HD 36526 




1.9 e-1 


3.0 


> 4.1 c6 


> 45 


< 3.9 


< 6.4 


> 1.18 









T Notes in Appendix. 
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hot to cool side of the so-called "bi-stability" jump at 
T efl « 25, 000 K (see below). 

For mass-loss rates, we also use the recipe given by 
Vink et al. (2000, 2001), assuming solar metallicity for all 
stars. This predicts an associated strong mass-loss increase 
of nearly an order of magnitude from the hot to cool side 
of this bi-stability jump, because iron recombination makes 
available more efficient driving-lines and so produces an in- 
crease in the line force. But note that, whereas the expected 
decrease in Voo over this bi-stability jump is empirically 
quite well established, this predicted increase in Mb=o is not 
yet observationally confirmed (e.g. Markova & Puis 2008). 

For comparison, we therefore also compute mass-loss 
rates based on the standard (finite-disk-corrected) Castor 
et al. (1975, hearafter CAK) scaling. Using the notation from 
Gayley (1995), this can be written in the form, 



Alt 



Q'cff 



L* ( QV e 



(1 + ttci) 1 /""" 1 - a cS c 2 



-l + l/a eft 



(7) 

where we adopt Q = 1000 and a e g ~ 0.55 for the full sample, 
to represent the normalisation and effective power-exponent 
of the line opacity distribution, where the latter has been 
adjusted to account for ionisation effects (Puis et al. 2000), 
and is in good agreement with the observationally inferred 
value for non-magnetic O-type stars (Repolust et al. 2004). 

The left panel of Figure 4 compares the two mass-loss 
rate values for our full sample of magnetic massive stars, 
while the right panel illustrates the shift in the confinement- 
rotation diagram resulting from switching between the two 
scalings. 

For the hotter O-type stars the CAK scaling agrees 
quite well with the Vink et al. recipe, and, because of the 
weak Ra ~ M^lg 4 dependence, this translates to a negligi- 
ble shift in the confinement-rotation diagram. Of course, this 
comparison only reflects uncertainties due to different theo- 
retical mass-loss descriptions. But recent multi-wavelength 
spectroscopic studies aiming to derive mass-loss rates in the 
O-type star domain typically yield rates that deviate from 
the Vink et al. prescription only by factors of ~2-3, if small- 
scale wind inhomogeneities ( "clumping" ) are adequately ac- 
counted for (Sundqvist et al. 2011; Najarro et al. 2011; 
Bouret et al. 2012) And as illustrated by Figure 4, such dis- 
crepancies barely affect stellar positions in the confinement- 
rotation diagram. 

For B-type stars, however, mass-loss rate differences are 
generally much larger, an order or magnitude or more near 
the bi-stability jump. Empirical mass-loss determinations 
for B-type dwarfs (which comprise most of our magnetic 
sample) are difficult at most, but studies of B-type super- 
giants have found a decrease in wind momentum compared 
to the theoretical Vink et al. predictions for the complete 
low-temperature region (Markova & Puis 2008). Further de- 
viations from theoretical wind momentum of similar mag- 
nitude have also been observed for some late O-type stars 
with so-called weak winds (for an overview see Puis et al. 
2008). 

Even with the weak M^Lq scaling, the shift in Ra as- 
sociated with these large deviations can approach 0.3 dex. 
For other quantities, such as the stellar spindown time, 
which scales as i? A ~ i/Vm b=0 , there can be a substantial 
change, by a factor of several, for different mass-loss values 



near and below the bi-stability region, as discussed further 
in §3.4. 

In summary, these relatively large systematic differences 
in the adopted mass-loss rate will be a important concern for 
performing detailed modelling of magnetosphere signatures 
for individual stars. However, it can be seen from Figure 4 
(right) that despite these large differences, the overall ap- 
pearance of the rotation-confinement diagram is not much 
affected and the basic, qualitative classification results pre- 
sented here are quite robust against errors in the wind pa- 
rameters. To maintain a uniform standard, all the presented 
magnetosphere parameter values in Table 6 are based on the 
Vink et al. scalings. 



3.3.2 Rotational oblateness 

As mentioned in § 3.1, calculation of the Alfven and Kepler 
radii requires the actual equatorial radius of the star, in prin- 
ciple accounting for any rotationally induced oblateness. In 
practice, a 15 percent oblateness requires Q. ~ 0.8, equivalent 
to W = 0.6, and so except for the most rapid rotators, the 
difference between the polar and equatorial radii is generally 
much smaller than the uncertainty in the radius determina- 
tion. In our sample, only three stars have a period short 
enough for oblateness to become potentially significant. For 
HD 182180 (ID 45) and HD 142184 (ID 47), we use the equa- 
torial radii derived by Rivinius et al. (2012) and Grunhut 
et al. (2012a) from spectral analysis including the oblate- 
ness. As no such analysis is available for HD 35502 (ID 57), 
we use the radius derived from the Stefan-Boltzmann equa- 
tion. 

For simplicity, we ignore the effect of gravity darkening 
on the wind driving from the stellar surface. For aligned ro- 
tators, the wind feeding the equatorial magnetosphere origi- 
nates from mid-latitudes, where gravity darkening is weaker. 
For non-aligned rotators, the magnetosphere will have a 
complex 3D structure that requires detailed modelling for 
each case. But in general terms, the maximum density occurs 
near Rk along the line defined by the intersection between 
the magnetic and rotational equators (Townsend & Owocki 
2005; Townsend et al. 2005). In this context, the relative con- 
finement and centrifugal support of the such magnetospheres 
should be well characterised by the Alfven and Kepler ra- 
dius relative to the star's equatorial radius, accounting for 
any rotational oblateness. 



3.3.3 Uncertainties in Kepler and Alfven radii 

Let us now explore the effect of stellar parameter uncertainty 
on the position of the stars in the confinement-rotation di- 
agram. As the radius and mass of the stars are generally 
derived from T e g, log(L*/L0) and log</, we propagate the 
uncertainty on these quantities, assuming they are indepen- 
dent. 

The Kepler radius can be expressed as, 



K K 9 -teg 
R* Ll/6 W 2/3' 



(8) 



In general, the rotational periods are accurate at 1-2 per- 
cent, so their uncertainty can be neglected. The quantities 
T c ft, log<j and log(L*/L0) have typical uncertainties of 10 




Figure 4. Left: Comparison between the mass-loss rate calculations by Vink et al. (2000, 2001) and the CAK scaling law described in 
§ 3.3.1. Note the large change in mass-loss over the bi-stability jump in the Vink et al. rates. Right: Shift in Alven radius from switching 
between the Vink et al. rates (filled symbols) and the CAK scaling law (empty symbols). The error bar in the lower right represents 
uncertainty in i?x and Ra estimated from the propagation of typical uncertainty in the stellar parameters only (~ 25 percent). The 
uncertainty in Ra is in fact dominated by systematics in the mass-loss rate determinations. However, the relative position in the magnetic 
confinement-rotation diagram, and hence the magnetospheric classification, is not too sensitive to these systematics, as described in § 3.3.1. 



percent, 0.2 dex and 0.25 dex, respectively. From equation 
8, they contribute 15, 15 and 10 percent uncertainties to 
Rk, for a total uncertainty of 23 percent. Figure 5 shows 
a histogram of the Kepler radius uncertainty distribution 
(grey shade) for our sample, confirming that the mean uncer- 
tainty is around 25 percent. A corresponding vertical error 
bar is shown in the confinement-rotation diagram in Figure 
4 (right). 

For the error propagation in Alfven radius, we follow 
the CAK scaling (for fixed a e ff = 0.55), for the dependence 
of wind momentum on stellar parameters (equations 6 and 
7), 



Ms=oVtx 



T- 



(9) 



With the typical uncertainties quoted above, T e ff, log<? and 
log(L*/L0) contribute respectively 22, 14 and 72 percent 
uncertainties to the wind momentum, for a total uncertainty 
of 76 percent. This uncertainty from stellar parameters is 
much smaller than that associated with the systematics dis- 
cussed in § 3.3.1. As such, we estimate the total uncertainty 
in Ra through the scaling, 



1/4 



Ra 



R* T eff (MjfcoKx,) 1 / 4 



(10) 



As a large fraction of our sample has only lower limits on 
the dipole field strength, let us ignore for now its contribu- 
tion to the uncertainty. Again with the uncertainties quoted 
above, T a s and log(L*/L ) contribute to 10 and 14 percent 
uncertainties, whereas Mb=oVoo contributes 20 percent un- 
certainty, for a total uncertainty in _Ra of 26 percent. Fig- 




20 40 
Radius uncertainty (%) 



60 



Figure 5. Distribution of uncertainties in our determination of 
the Kepler radius (shaded) and Alven radius (hatched), estimated 
by propagating the uncertainty on the stellar parameters (T e g, 
log(L*/L ) and logs). 



ure 5 presents the uncertainty distribution attributed to the 
stellar parameters for Ra (hashed histogram) with a mean 
value again confirming the above estimate. One can see that 
the corresponding horizontal error bar in the confinement- 
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rotation diagram of Figure 4 (right) is smaller than the sys- 
tematic uncertainty associated with the theoretical mass- 
loss rate determination. 

Moreover, since the dipolar field strength is generally 
constrained with an accuracy of ~ 30 percent, it would 
only contribute a 15 percent uncertainty to Ra, again much 
smaller than the systematic uncertainty from the wind mo- 
mentum. 



3.4 Spindown time and spindown age 

Let us next turn to considering the rotational evolution for 
our sample of magnetic massive stars. The angular momen- 
tum loss rate from a magnetised wind can be written in 
terms of the mass-loss rate, the Alfven radius Ra, and the 
stellar rotation frequency to = V TO t/R* (Weber & Davis 
1967; ud-Doula et al. 2009), 



(11) 



The associated timescale for magnetic-wind-induced spin- 
down of the stellar angular momentum J = Iui can then be 
written in the form 



T.T 



J 
j 



3 fM,R 2 ,u 
2 M b=0 RIcj 



= n f™ 



TJB=0 



R*^ 
Ra 

R*_ 
Ra 



(12) 



where tm = M*/Mb=o is a characteristic mass-loss 
timescale, and tj,b=o defines the spindown time in the case 
of no magnetic field (i.e. Ra — R*)- The star's moment of 
inertia I = /M*.R* can be evaluated from the radius of gy- 
ration j3 — f 1 / 2 tabulated from internal structure models 
such as Claret (2004). If we assume for simplicity a fixed 
radius R, and moment of inertia factor / « 0.1, as well as 
a constant angular momentum loss rate J, then the stellar 
rotational period P will simply increase exponentially with 
age t from its initial value P a , 

P e t/T -> . 



P(t) 



(13) 



We can then use equation 13 to define a star's spindown 
age, t s , in terms of the spindown time rj, and its inferred 
present-day critical rotation fraction W = P rh/P relative 
to its initial rotation fraction W at age t = 0, 

' s = lnW„-lnW. (14) 



Taking the initial rotation to be critical, W 
simple upper limit to the spindown age, 



ln(l/VK) . 



1, yields a 



(15) 



If the initial rotation is subcritical, W < 1, then the actual 
spindown age is shorter by a time At s = rj lnW^o. 

As noted previously, the extra axes in Figure 3 give the 
spindown timescale rj normalised by the value in a non- 
magnetised wind (i.e. by how much the magnetic braking 
enhances the stellar spindown) along the top, and the maxi- 
mum spindown age f s , max normalised by the spindown time 
(i.e. the number of spindown e-folds) along the right. For 
each of the individual magnetic OB stars, columns (8) and 
(9) of Table 6 also list estimated values for respectively the 
spindown time rj and the maximum spindown age i s ,max in 



Myr. Future studies can thereby compare t s ,max with other 
indicators of stellar age, for example from stellar evolution 
tracks or cluster association. To the extent that such inde- 
pendent age estimates are available, then within the limits of 
the stated assumptions of constancy in R*, f, and J, a com- 
parison with this spindown age could be used to estimate an 
initial rotation fraction W . 

More immediately, note that among the full magnetic 
sample, many of the most slowly rotating stars are O-type 
stars. The high luminosities of these stars drive strong stellar 
winds that lead to a rapid angular momentum mass loss and 
thus very short spindown times. These characteristics help 
to explain their very slow rotation relative to many of the 
B-type targets. Except for Plaskett's star (ID 6), which has 
likely been spun up by binary interaction and show CM- 
type emission at high velocity (Grunhut et al. 2012b), all 
the rapidly rotating stars near the top of Figure 3 are lower 
luminosity B-type stars with weaker winds; for magnetic B- 
type stars, the spindown time is thus generally longer than 
for the magnetic O-type stars, typically several Myrs. 

Indeed, extended photometric monitoring of the 
strongly magnetic B-type star a Ori E (ID 31) has provided a 
direct measurement of the change in rotation period, yield- 
ing a spindown time of 1.34 Myr (Townsend et al. 2010, 
2012). This is remarkably close to the spindown time of 
1.4 Myr predicted previously by the scaling developed from 
MHD simulations (ud-Doula et al. 2009), but such very close 
agreement was likely fortuitous given the uncertainties in the 
mass-loss rate and stellar parameters. Indeed, the Vink et al. 
(2000, 2001) mass-loss rate we use here is roughly a factor 10 
smaller than the CAK mass-loss rate assumed by ud-Doula 
et al. (2009), leading to a factor ~ \/W longer estimate for 
the spindown time, 4.6 Myr. This emphasises that our listed 
values spindown time and age are only estimates accurate 
to within a factor of 3 or so 7 . 



4 Ha AS A MAGNETOSPHERIC PROXY 

We now explore how magnetospheric Ha emission char- 
acteristics correlate with their position in the magnetic 
confinement-rotation diagram. 



4.1 Identification of Ha magnetospheric emission 

In Table 6, column (11) indicates the emission (em) versus 
absorption (abs) nature of the Ha line, and, if enough ob- 
servations (> 5) are available, whether the profile is stable 
(stab) or variable (var). We flag Herbig Be stars (HeBe) be- 
cause of the difficulty in disentangling magnetospheric emis- 
sion from the emission produced by the accretion disks char- 
acteristic to this class of pre-main sequence stars. Similarly 
in the case of spectroscopic binaries, slowly pulsating B-stars 
and /3Cep stars (indicated in column 3), variation in the 



7 A change of period has also been measured for HD 37776 (ID 34, 
Mikulasek et al. 2011) with a spindown timescale of 0.37 Myr, 
comparable to our rj of 0.68 Myr. However, as discussed by ud- 
Doula et al. (2009), the complex field geometry of HD 37776 
(Kochukhov et al. 2011) will have a potentially strong impact 
on the angular momentum evolution of this star. 
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absorption line profile may have non-magnetic origins, but 
emission can generally be attributed to a magnetosphere. 

This emission can have various distinguishing charac- 
teristics: (1) a central absorption core with broad emission 
wings that extend well beyond the photospheric v sin i (Fig. 
6, left), (2) strong, narrow emission that overwhelms the 
photospheric absorption profile (Fig. 6, right, grey profile), 
(3) weak overlying emission that only partially fills the un- 
derlying absorption (Fig. 6, right, black profile). 

For type (3), there can be confusion with the line- filling 
effect of a non-magnetised stellar wind; a clear identification 
requires monitoring for rotational modulation 8 . For type (1) 
and (2), a single observation can suffice to identify a mag- 
netospheric origin. 

For type (1) the extended emission wings suggest 
plasma held in extended rigid-body rotation around the star, 
presumably by the stellar magnetic field. These correspond 
to centrifugal magnetospheres, as described by Townsend & 
Owocki (2005). Multiple occurrences of this type of emis- 
sion can be found in Bohlender & Monin (2011), Oksala 
et al. (2012), Grunhut et al. (2012b), as well as the refer- 
ences listed in Table 2. 

For type (2), the narrow central emission suggests that 
the trapped plasma is kept at low velocities, without much 
broadening from rotation or from a high-speed outflow like 
in a non-magnetic stellar wind (Sundqvist et al. 2012). These 
correspond to dynamical magnetospheres. Example of such 
emission can be found in Howarth et al. (2007), Grunhut 
et al. (2012c), Wade et al. (2012a), as well as the references 
listed in Table 2. 



4.2 Ha emission in the classification diagram 

Figure 7 (left) again plots stars in the confinement-rotation 
diagram, with symbols now coloured to mark the presence 
(dark pink) or absence (light purple) of magnetospheric Ha 
emission. Herbig stars are omitted here because of their in- 
trinsic emission not associated with magnetic fields. While 
stars with and without emission are found throughout the 
diagram, note that in the DM region with _Ra < Rk, all 
the emission occurs, (with just one exception, HD 156424; 
ID 36) in O-type stars, for which the large luminosity sug- 
gests the wind feeding of the DM is strong enough to build 
up sufficient density for emission within the dynamical in- 
fall timescale. The same strong wind mass loss that feeds the 
DM emission means that they have relatively strong angular 
momentum loss that spins down the stars to their observed 
slow rotation rates near the bottom of the confinement- 
rotation diagram. 

Conversely, in the CM region with Ra > Rk, the emis- 
sion (again with one exception, Plaskett's star, ID 6) occurs 
in B-type stars, for which the lower luminosity and wind 
feeding requires the longer retention timescale of a CM to 
build up sufficient density for emission. In fact, all the non- 
emitting stars in this region are also B-type, indicating that 



8 Modulated variations in the core of absorption lines could also 
have other origins, for example changes in the photospheric struc- 
ture due to large helium abundance inhomogeities on the surface 
of chemically peculiar stars (e.g. a Cen; ID 63, Bohlender et al. 
2010) 



a CM is a necessary but not sufficient condition for emis- 
sion for such low luminosity stars with relatively weak wind 
mass loss. Most of the B-type stars with emission are in the 
extreme upper right of the diagram, with both strong con- 
finement and rapid rotation. Their wide separation from the 
Ra = Rk line implies a large radial extent for their CM. 

Overall, this link between Ha emission and location in 
the magnetic confinement-rotation diagram provides a use- 
ful categorisation that connects the rotation, mass loss, and 
circumstellar emission properties of massive star magneto- 
spheres. 

4.3 Magnetospheric vs stellar properties 

To explore further this categorisation, Figure 8 plots (again 
for all the non-Herbig stars) the log of the ratio Ra/Rk 
(column 10 of Table 6) vs. stellar effective temperature 
T e ff (left) or bolometric luminosity (right), with symbols 
again marking spectral type, coloured for the presence (dark 
pink) or absence (light purple) of magnetospheric Ha emis- 
sion (Figure 9 provides a finding chart). The single upward 
arrows indicate stars that could be shifted upward due to ei- 
ther a higher polar field than the minimum inferred from the 
available longitudinal field measurements (increasing -Ra), 
or a higher rotation rate than the minimum inferred from 
the measured vsini (decreasing Rk)', the double upward ar- 
rows indicate stars for which both limits are at play. 

The solid horizontal line at Ra = Rk separates the do- 
mains for dynamical magnetospheres at the bottom from 
centrifugal magnetospheres at the top, with the distance 
above the line characterising the radial extent for the cen- 
trifugal support. Each plot again shows that Ha emission 
occurs both in O-type stars to the left, and in B-type stars 
to the (mostly upper) right; but the demarcation is partic- 
ularly distinct in the plot vs. bolometric luminosity. 

In that plot, the vertical dashed line corresponds 
roughly to the main sequence transition from O- to B-type 
stars (Martins et al. 2005). The O-type stars to the left all 
have clear Balmer emission. Except for Plaskett's star (ID 6), 
which has likely been spun up by binary interaction, they 
also are all relatively slow rotators with Ra < Rk (DM). 
Their high luminosity means they have strong stellar winds, 
implying both a rapid stellar spindown, as well as a suf- 
ficient magnetospheric density to give line emission in the 
short residence time for a dynamical magnetosphere. 

For B-type stars to the right of this vertical line, emis- 
sion is most common in the most rapidly rotating stars above 
the horizontal dotted line at log(i?A/RK) = 0.9. The 4 stars 
above this line without detected emission are all relatively 
late type stars, with low luminosity and so likely a very low 
wind mass-loss rate to feed the expected CM. The 3 stars 
below this line with detected emission have arrows indicat- 
ing they could shift upward, with the two lower luminosity 
stars (HD 156424 and HD 156324; ID 36, 37) having dou- 
ble arrows indicating potentially significant revision in both 
field strength and rotation rate. Indeed, although the cur- 
rent MiMeS observations do not allow for clear determina- 
tion of a rotation period, in both cases the nightly variation 
of longitudinal field measurements points toward periods of 
order a day. These stars are prime candidates for follow up 
observations. 

The third, relatively high-luminosity B-type star 
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Figure 6. Example of Ha line profiles, from MiMeS observations, for a centrifugal magnetosphere (left, a Ori E, ID 31) and a dynamical 
magnctosphere (right, HD 191612, ID 4) and minimum (black) and maximum (grey) emission. Note how the extended emission wings 
are located at large velocities (outside usini) for the CM, whereas the central emission is localised in a narrow range of velocities (inside 
v sin i) in the line core for the DM. 




Figure 7. Location of magnetic massive stars in the magnetic confinement-rotation diagram (see Figure 3). The symbols are coloured 
to mark the presence (dark pink) or absence (light purple) of a magnetospheric signature in Ha emission (left) and UV resonance line 
(right), as described in Table 6, and in §4.1 and §5.1 respectively. The symbols are empty when no information is available. 



CMa; ID 15) has a well determined period and Kepler 
radius, but still has a single arrow from the limited polar 
field estimate. Its current position - just above the horizon- 
tal solid line, and just to the right of the vertical dashed line 
- makes it a particularly interesting test case for magneto- 
spheric models, very near the transition from DM to CM, 
and from O-type stars to B-type stars mass loss. 



Indeed, if the CM/B-type star occurrence of Ha emis- 
sion depends on a combination of the radial extent of the 
CM (set by Ra/Rk and thus the vertical position in Figure 
8) and on the mass-loss rate feeding the CM (set by the lu- 
minosity and thus the horizontal position in Figure 8), then 
we can identify a possible further division along the illustra- 
tive diagonal dot-dashed line in Figure 8. It would thus be 
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Figure 8. Location of magnetic massive stars in a log-log plot of Ra/^k vs - the effective temperature (left) and the luminosity (right). 
The symbols are coloured to mark the presence (dark pink) or absence (light purple) of magnctosphcric Ha emission, as described in 
Table 6 and §4.1. The symbols are empty when no Ha information is available. Single arrows indicate a limit on either R\ or i?Ki 
whereas double arrows mark stars for which both and i?x are limits. In the righthand diagram, the vertical dashed line represents 
the luminosity transition between O-type and B-type main sequence stars. The horizontal dotted line and the diagonal dot-dashed line 
are illustrative division of the CM domain according to potential mass leakage mechanisms (see discussion in §4.3). 



< 



ao 
o 



□12 



* 
^9 



□11 



45. A7^ 7 

,,3430*49 otf c , 
«62 

22*^ „*46 
,24 50 4 » 42 

^20 ^35 "aa 

□14 
□13 



□19 £7 



4.6 4.5 4.4 4.3 

log( r e „ ) (kK) 



4.2 



o 



#57 .45.47 

"26 31 ^Igl 1 55 64; 

34^59 58 6 o°- 
30 .62 
»28 22 .51 

"25 .46 
.42 24ji 5o 




log (L./L ) 



Figure 9. Finding charts for the location of magnetic massive stars in a log-log plot of /?a/^?k vs. the effective temperature (left) and 
the luminosity (right). The label numbers correspond to the ID in column (1) of Table 1. 



of particular interest to clarify the position, and the emis- 
sion properties, of stars with current placements near this 
illustrative diagonal line. 

Establishing empirically whether the onset of emission 
is better delineated by the horizontal dotted line or the 
diagonal dot-dashed line has potentially important impli- 



cations for our theoretical understanding of the magneto- 
spheric mass budget. The former would indicate that the 
CM mass depends mainly on the capacity for the magnetic 
field to confine centrifugally supported material, which even- 
tually fills to a fixed level even if at the slow rate from a weak 
wind. The latter would indicate a competing leakage from 
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the CM that decreases with distance above the _Ra = Rk 
line. To build up sufficient density for emission, stars near 
the Ra = -Rk line with high leakage require a high mass-loss 
rate and thus high luminosity, representing the left end of 
the diagonal; stars further above the line with lower leakage 
can fill their CM even with the weaker wind from a lower 
luminosity, representing the upper right end of the diagonal. 

As a shorthand, we might identify these as the 'capacity' 
vs. 'leakage' models for determining the onset of CM emis- 
sion. Hopefully, the classification and physical arguments 
here will motivate a concentrated observational program to 
clarify the position, and the emission properties, of the key 
stars for establishing this discrimination. 



5 OTHER MAGNETOSPHERIC PROXIES 
5.1 Ultraviolet variability 

In hot, massive stars, strong UV resonance lines like 
C iv AA 1548, 1550, Si iv AA 1393, 1403 and N v AA 1239, 1243 
are typically used as diagnostics of the stellar wind. In O- 
type stars with dense winds, the line profiles generally ex- 
hibit a characteristic P-Cygni profile showing red-side emis- 
sion and blue-side absorption, with the blue edge of the lat- 
ter marking the wind terminal speed. In B-type stars with 
weaker winds, the emission is weak or absent, and the blue 
edge of the shallower absorption may not extend to the 
terminal speed. Both types can exhibit intrinsic variability, 
but this is most common, distinctive and well-studied in O- 
type stars, for which it is generally characterised by discrete 
absorption components (DACs) that start near line-centre 
and slowly propagate across the blue absorption trough (e.g. 
Howarth & Prinja 1989; Kaper et al. 1996). These are likely 
representations of spiral-shaped density compressions, re- 
ferred to as Corotating Interacting Regions (CIRs, Mullan 
1986), caused by faster moving streams overtaking slower 
moving streams, where the interacting interface between the 
two is shocked. The projected velocity in the line of sight 
progresses because of the stellar rotation, rather than be- 
cause of the outflow itself (Cranmer & Owocki 1996) 

In magnetic OB stars these UV lines can be strongly af- 
fected by changes in velocity, density, and/or ionisation bal- 
ance. Indeed, it was recognised early on that periodic varia- 
tions of the resonance lines could point toward the presence 
of a rotating magnetosphere (e.g. Shore & Brown 1990; Hen- 
richs et al. 1993; Walborn & Nichols 1994). In contrast, the 
appearance of DACs is often found to be cyclical but never 
strictly periodic. Moreover, unlike the blueward-propagating 
DACs, in magnetic OB stars UV line variation occurs nearly 
coherently and synchronously over the full velocity range of 
the profile (e.g. Marcolino et al. 2012; Henrichs et al. 2012). 
Thus, even in those stars without sufficient monitoring to 
clearly establish a period, one could use the morphological 
character of variations between two or more observations to 
flag the likely presence of a strong field (e.g. Henrichs et al. 
2012). 

In B-type stars, UV profiles show only shallow (if any) 
blueward wind absorption and weak or absent redward emis- 
sion. The appearance of strong redward emission along with 
filling in of the absorption (see Figure 5 of Shore & Brown 
1990) can likewise be used to flag the likely presence of a 



magnetosphere, even without multiple observations to show 
variability. 

Column (12) of Table 6 gives a summary characterisa- 
tion of UV resonance lines for the full sample of magnetic 
stars. The listed UV signatures of a field include periodicity 
(per), profile variability with morphology similar to periodic 
stars (var), and, for B-type stars, distinct redward emission 
with missing blueward absorption (em) . Stars lacking a clear 
UV magnetic signature are those with only pure absorption 
(abs), and those with 5 or more observations showing sta- 
ble absorption (stab abs). When available, these character- 
isations are from the literature, as summarised in Table 2, 
and otherwise are based on visual inspection of IUE archive 
spectra. 

In analogy with the organisation of Ha signatures 
shown in the left panel of Figure 7, the right panel again 
plots stars in the magnetic confinement-rotation diagram, 
but now with stars showing one or more UV signatures for 
a field marked in dark pink, and those with absorption pro- 
files consistent with lack of a field marked in light purple. 
Stars without UV observations (or IUE spectra with too low 
signal-to-noise ratio) have empty symbols. 

Note that, in contrast to Ha emission, such UV mag- 
netic signatures occur throughout the diagram, and for all 
spectral types. In particular, B-type stars with weak winds 
show a UV magnetic signature even in the slow-rotation, 
DM region, for which B-type star Ha emission is not seen. 
Thus UV variation seems to be a wide-spread phenomenon 
among magnetic OB stars, as long as some confinement is 
present, and therefore represents a relatively robust proxy 
of magnetism. In fact, a number of magnetic OB stars had 
been first identified as peculiar UV stars (e.g. Henrichs et al. 
1993; Neiner et al. 2003; Henrichs et al. 2012). 

The few stars without signs of UV variability cluster 
at lower 7?a, but many other stars show UV field signa- 
tures in the same region of the diagram. As suggested by 
Shore & Brown (1990) for some of the Bp stars, the exact 
behaviour of the variability might be closely tied with the 
geometry of the magnetic field with respect to the observer. 
Thus detailed modelling of the UV line profiles for magnetic 
OB stars may help constrain the geometry of the magneto- 
spheres and clarify the velocity and ionisation structure of 
the trapped material. 

Indeed, UV resonance-line synthesis from MHD mod- 
els shows clear P-Cygni absorption troughs that are mod- 
ulated on the rotation phase. For relatively strong lines, 
such troughs are actually deeper when viewing the mag- 
netosphere pole-on than equator-on (ud-Doula 2008). This 
somewhat counterintuitive effect occurs because the over- 
dense material around the magnetic equator is characterised 
by very low velocities, whereas the outflow above the pole 
more closely resembles that of a normal, non-magnetic wind. 
Thus the absorption column above the pole covers a much 
wider velocity range, leading to wider and deeper troughs. 

However, further calculations also suggest that the 
phase variability of such UV lines is quite sensitive to the ac- 
tual strength of the line itself, and so may depend critically 
on the stellar mass-loss rate as well as on the magnetospheric 
ionisation state (Sundqvist et al. in prep). The strong UV 
lines in HD 108 (ID 7) indeed seem to display the charac- 
teristic variability described above (Marcolino et al. 2012), 
but those in 9 1 Ori C (ID 3) show effectively the opposite 
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behaviour (ud-Doula 2008). Thus further modelling work is 
still needed to fully understand how magnetic fields affect 
the formation of UV lines of OB stars. 

5.2 X-rays 

Massive stars are generally X-ray bright due to the intrin- 
sic instability of the line-driving mechanism for radiative 
stellar winds (Feldmeier et al. 1997; Runacres & Owocki 
2002; Dessart & Owocki 2003), with a well known canonical 
value for early OB star X-ray luminosity, Lx, of ~ 10 _7 ' 2 L* 
(Berghoefer et al. 1997; Naze et al. 2011; Gagne et al. 2011). 
The magnetically channeled wind shocks (MCWS scenario) 
associated with magnetic massive stars should also gener- 
ate even stronger and harder X-ray emission, by the radia- 
tive cooling of the shock heated plasma in the magneto- 
sphere (Babel & Montmerle 1997a). For example, the X-rays 
from the O-type star 9 1 Ori C (ID 3) are more luminous and 
harder than in typical O-type stars, and modulated by the 
rotational period. Gagne et al. (2005) used 2D MHD simu- 
lations, including an explicit energy equation, to track the 
shock heated material and its radiative cooling, and were 
able to reproduce the X-ray properties of 6 1 OriC, includ- 
ing the star's elevated X-ray luminosity, high plasma tem- 
perature, rotational modulation, and narrow spectral lines. 
Therefore, it seems at first glance that luminous, hard and 
variable X-ray emission could be a proxy for magnetism in 
massive stars. 

However, these characteristics are not always present in 
magnetic massive stars. For example, the B-type star r Sco 
(ID 11) is X-ray luminous and indeed displays a hard X-ray 
spectrum (Mewe et al. 2003), but it does not show evidence 
of rotational modulation (Ignace et al. 2010). The B-type 
star NU Ori (ID 12) does not show any significant variability 
over the duration of a ~ 10 d Chandra observation (Stelzer 
et al. 2005), and has a soft spectrum. Another prototypi- 
cal magnetic O-type star is the Of?p star HD 191612 (ID 4) 
which is quite luminous, but has a rather soft spectrum 
(Naze et al. 2010). Oskinova et al. (2011) recently exam- 
ined a small subset of magnetic B stars and noted that they 
too have diverse X-ray properties, including a few that are 
not X-ray overluminous at all. 

We present here a first attempt to cast the X-ray charac- 
teristics of our large sample of magnetic OB stars as a func- 
tion of their magnetospheric properties, focusing on just the 
X-ray luminosity. We perform a review of the literature to 
extract X-ray fluxes for the stars in our sample. Where pos- 
sible, we use X-ray fluxes derived from pointed observations 
by modern X-ray observatories (Chandra and XMM) and 
reported in papers that carefully model the emission prop- 
erties, correcting for interstellar absorption. A large majority 
of the O-type and very early B-type stars in our sample fall 
into this category. And for these stars, differences in the in- 
strument bandpasses and uncertainties associated with the 
multi-temperature emission modelling and the ISM correc- 
tion should lead to errors in the reported X-ray fluxes of 
less than a factor of two. We correct all of the literature X- 
ray luminosities for the distances adopted by the authors of 
each paper to derive an X-ray flux, and then recompute the 
X-ray luminosity using the distances we adopt for each star, 
which of course are consistent with the distances we use for 
the bolometric luminosity determinations. We then compute 



the X-ray efficiency ratio log(Lx/£*) (column 13 of Tabic 
6), so that even if better distance determinations are made 
for some of these objects in the future, their log(Lx/-I/*) 
values will not have to be adjusted. 

For many of the later B-type stars, no X-ray measure- 
ments exist, and for others only survey data, primarily from 
ROSAT, exists. The X-ray fluxes derived for these stars are 
more uncertain, primarily because of the lack of detailed 
spectral modelling and in some cases the lack of detailed 
ISM absorption corrections. Additionally, the bandpass of 
ROSAT is softer than that of either XMM or Chandra, fur- 
ther skewing comparisons between the derived X-ray lumi- 
nosities. There are similar considerations for the small num- 
ber of B-type stars for which only EINSTEIN measurements 
exist (Grillo et al. 1992). A more conservative estimate of the 
log(Lx/£*) uncertainties for these stars is required, with the 
overall error being probably up to 0.5 dex. Another factor 
potentially affecting our X-ray luminosity determinations is 
the contribution from unresolved binary companions (e.g. 
Petit et al. 2012). This phenomenon is more likely to be 
important for later B-type stars, observed with X-ray tele- 
scopes with poorer spatial resolution, and with lower intrin- 
sic X-ray luminosities such that low-mass PMS companions 
could account for much of the observed X-ray emission for 
a given star. However it is unlikely that all the X-ray bright 
magnetic B-type stars are affected by binarity. For example, 
Gagne et al. (2011) have shown that the pre-main sequence 
population of the Carina Complex cannot explain all the X- 
ray emission of B-type stars and that some of them must be 
intrinsically X-ray bright. 

In Figure 10 (left) , we plot the stars in the Ra-Rk plane 
with a colour coding representing the X-ray efficiency ra- 
tio in bins of 0.5dex ([>-6.5], [-6.5,-7.0], [-7.0,-7.5] and 
[< — 7.5]). The dark pink shades are for stars with X-ray ef- 
ficiency greater than the canonical value of log(Z/x/£*) 

7,0 for O-type stars. All the O-type stars show some level 
of overluminosity (log(Lx/i*) > —6.7). Some of the B-type 
stars also show overluminosity. Most of them are located 
in the upper part of the CM region, although a few over- 
luminous, very early B-type stars are located in the DM 
region. The right panel of Figure 10 presents a logarithmic 
plot of Ra/Rk vs. the bolometric luminosity. One can see 
that enhanced X-ray emission generally occurs for the most 
luminous magnetic OB stars. The B-type stars with inter- 
mediate luminosities seem to have X-ray emission typical 
for their spectral type, although it has been shown that the 
Lx-L* relation breaks down at spectral type B2 and that 
the typical Lx is W~ s L ir or lower for later spectral types 
(Cohen et al. 1997), implying that any mid B-type star in 
one of the three highest log(Lx/i*) bins is overluminous for 
its spectral subtype. Finally, some low luminosity stars with 
large Ra/Rk (extended centrifugal magnetospheres) show 
enhanced X-ray emission. 

A potential explanation for the X-ray emission enhance- 
ment in CMs is the centrifugal acceleration for fast rotators, 
which contributes by propelling the confined material up the 
magnetic loop leading to stronger shocks than what could be 
achieved by radiative acceleration alone like in a DM. How- 
ever, Rigid-Field Hydrodynamical (RFHD; Townsend et al. 
2007) simulations predict that the overall X-ray flux of stars 
in the CM region is also quite sensitive to the mass-loss rate. 
The distribution in temperature of the differential emission 
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Figure 10. Location of magnetic massive stars in the magnetic confinement-rotation diagram (left) and in a log-log plot of Ra/Rk 
vs. the luminosity (right). The stars are colour-coded according to their X-ray luminosity with respect to their bolometric luminosity 
(log(Lx/^*))i i n bins of 0.5dex. The dark pink shades represent stars with log(Lx/^*) greater than the canonical value of — 7 for early 
OB stars and the light purple shades for stars below this value. The symbols are empty when no X-ray information is available. 



measure (DEM) is governed by both the pre-shock and post- 
shock characteristics of the magnetosphere. Both of these are 
affected by the wind properties, with the post-shock cooling 
length being longer for lower-density wind flows, leading to 
softer and weaker emission (Hill et al. 2011). 

Therefore, a complete survey of the X-ray properties of 
the magnetic OB stars would be highly desirable as X-ray 
emission could provide a different perspective on the struc- 
ture and dynamics of magnetospheres, and the shock physics 
occurring in both DMs and CMs. Future studies should in- 
clude a consistent and uniform analysis of (1) X-ray plasma 
temperature distributions and (2) time- variability. Although 
a re-analysis of all the available X-ray observations is be- 
yond the scope of the current paper, the work presented 
here can be used as a starting point for identifying inter- 
esting stars for which X-ray data already exists, as well as 
identifying stars with interesting positions in the magnetic 
confinement-rotation diagram for which acquiring X-ray ob- 
servations should be a priority. Another key development 
for understanding the trends identified above in the X-ray 
emission of magnetic massive stars would be more accurate 
and secure wind mass-loss rates for the B-type stars. 



6 CONCLUSIONS 

The Magnetism in Massive Stars (MiMeS) project aims to 
study the scope and impact of stellar magnetism in mas- 
sive stars using high resolution and high signal-to-noise ra- 
tio spectropolarimetric observations from Large Program 
time allocations. Within this context, the present study 
had two main goals: (1) To compile an exhaustive and 
well-documented list of confirmed magnetic, hot OB stars 
that are directly detected through the Zeeman effect, and; 



(2) To organise the stars in a way that accounts for both 
the strength of magnetic confinement of the stellar wind 
(through ?7» or Ra) and the dynamical role of stellar rota- 
tion (through W or -Rk). Key results are: 

• We have provided a compilation of relevant stellar pa- 
rameters for our magnetic sample. We used the luminosity, 
mass, and radius obtained from modern spectral modelling 
from the literature, when available; otherwise, these were de- 
rived from a classical bolometric correction approach, and 
from SED fitting with the code Chorizos when a complete 
set of photometry was available. We have also compiled rota- 
tional periods and dipolar field strengths, as well as binarity 
and pulsation status. 

• Using these parameters, we have placed the full sam- 
ple of magnetic stars in a classification plane, the magnetic 
confinement-rotation diagram, characterising stellar rota- 
tion (as Rk or W) vs. wind magnetic confinement (as Ra 
or rj t ). 

• We identified key domains within the magnetic 
rotation-confinement diagram, representing weakly magne- 
tised winds with r]* < 1, or dynamical magnetospheres (DM, 
with J?* < Ra < Rk) vs. centrifugal magnetospheres (CM, 
with R t < i? K < Ra)- 

• We have associated Ha line emission characteristics 
with position in the confinement-rotation diagram. Slowly 
rotating O-type stars show DM magnetospheric emission, 
in contrast to B-type stars which generally only show CM 
magnetospheric emission. 

• In a plane plotting the ratio Rk/Ra vs. stellar lumi- 
nosity, we found a clear separation between O-type star DM 
emission and B-type stars for which appearance of CM emis- 
sion requires higher Rk / Ra for lower luminosity stars. This 
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suggests that the CM leakage mechanism may depend on 
the degree of magnetic confinement. 

• We have also associated other magnetospheric proxies 
with position in the confinement-rotation diagram. UV res- 
onance line variation occurs in all magnetosphere regimes 
and for stars of all temperatures; although detailed mod- 
elling will be needed in the future, UV spectroscopy seems 
an excellent proxy for identifying new magnetic OB stars. 
The earliest magnetic OB stars generally show X-ray over- 
luminosity, as do the low-luminosity B-type stars with large 
centrifugal magnetosphere volumes (high Ra/Rk)- 

• We have calculated magnetic spindown timescales (tj) 
and inferred spindown ages (t B ) for each star in our sam- 
ple. O-type stars with strong winds have short spindown 
timescales and so mostly are slow rotators located in the 
DM region; B-type stars with weaker winds have longer spin- 
down timescales, and thus extend well into the CM regime 
with rapid rotation. 

• Finally, we have identified stars which will be prime can- 
didates for follow-up studies (with either unknown periods 
or only dipole field strength lower limits) that would lead 
to a more accurate placement on the confinement-rotation 
diagram, hence providing more clues to the answers of some 
of the questions raised in this paper. 
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APPENDIX A: NOTES ON INDIVIDUAL 
STARS 

This section contains additional remarks about certain stars, 
marked by a dagger in Table 1, concerning our choice of 
parameters, providing alternative calculation in case of dis- 
agreements in the literature, or specific information about 
binarity and other relevant characteristics. The sections are 
numbered according to the ID number of each stars, as given 
in column (1) of Table 1. 

A3 HD 37022 (6 1 OriC) 

6 OriC is a single-lined spectroscopic and astrometric bi- 
nary with an 11-yr period and e « 0.6 (Kraus et al. 2009). 

A4 HD 191612 

HD 191612 is a double-lined spectroscopic binary with P — 
1548 d and e = 0.5, and a early B-type companion (Howarth 
et al. 2007; Wade et al. 2011b). 

A6 HD 47129 (Plaskett's star) 

Plaskett's star is a high-mass (M to tsini = 93M©) O+O 
double-lined spectroscopic binary with a short period (P — 
14 d) circular orbit (Linder et al. 2008). The magnetic field 
is associated with the rapidly rotating (vsini = 300km s _1 ) 
secondary star (Grunhut et al. 2012b). 

A8 ALS 15218 (TV 16-22) 

We use the effective temperature and luminosity derived 
from the cluster photometry analysis of Gagne et al. (2011). 
We assume log(g) = 4.0. 

A10 HD 37742 (COriAa) 

HD 37742 is an astrometric double-lined spectroscopic bi- 
nary (Hummel et al. 2000) with a 091b primary and an 
early B-type star companion. The preliminary dynamic mass 



derived by Rivinius et al. (2011) (2.48 ± 5.6 M ) is smaller 
than the mass derived by Bouret et al. (2008) (39 ± 8 M Q ). 

A comparison of the disentangled component spectra 
with the published Zeeman magnetic signature confirms that 
the signature cannot originate from component Ab, since its 
lines are too narrow (Th. Rivinius, priv. com). 

A12 HD 37061 (NU Ori) 

NU Ori is a double-lined spectroscopic binary with P = 19 d 
and e = 0.14 (Abt et al. 1991), with a magnetic early B-type 
primary and a late B-type companion (Petit et al. 2008). 

A18 HD 205021 (j3 Cep) 

P Cep is a double-lined spectroscopic binary. The magnetic 
primary (component A) is the prototype of a class of pul- 
sating hot stars. The secondary (component Aa) is an Ha- 
emitting classical Be star (Catanzaro 2008). 

A19 ALS 15211 (TV16-13) 

We use the effective temperature and luminosity derived 
from the cluster photometry analysis of Gagne et al. (2011). 
We assume \og(g) — 4.0. 

A20 HD 122451 (/? Cen) 

P Cen is a double- lined spectroscopic binary with P — 356 d 
and e = 0.8 (Ausseloos et al. 2006), and components of 
nearly identical mass. The magnetic field detection is as- 
sociated with the narrow-line primary star (Alecian et al. 
2011). 

A24 HD 96446 (V 430 Car) 

Matthews & Bohlender (1991) observed photometric vari- 
ations with a period of 0.85 d, as well as other shorter pe- 
riods interpreted as /3 Cep-type pulsations. However, other 
photometric periods are possible and compatible with the 
low-amplitude variations of the longitudinal field measure- 
ments (Neiner et al. 2012b). We use the long period of 5.73 d. 
The shortest period of 0.85 d would yield W = 0.46 and 
Rk = 1.7 R,. 



A30 HD 37017 (V 1046 Ori) 

V 1046 Ori is a double-lined spectroscopic binary with a 
18.6 d period (e = 0.4). The field detection is associated 
with the B2 He-strong primary (Bolton et al. 1998). The 
companion is a late B-type star. 



A32 HD 149277 

An inspection of archival HARPS and MiMeS ESPaDonS 
spectra revealed that HD 149277 is a double-lined spectro- 
scopic binary with the magnetic field detection associated 
with the lower vsmi component (Petit et al. in prep). 
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A34 HD 37776 (V 901 Ori) 

We choose the temperature and luminosity from Landstreet 
et al. (2007), in order to be consistent with the radius 
used by Kochukhov et al. (2011) in their magnetic analysis. 
Kochukhov et al. (2011) have shown that the field structure 
is much more complex than a dipole. Their surface magnetic 
field reconstruction displays a surface field varying from 5 kG 
to 30 kG across the stellar surface (see their Fig. 4). As their 
associated mean magnetic field modulus varies between 13 
and 16 kG, we use a dipolar strength of 15 kG to estimate 
the wind confinement. It is important to keep in mind that 
the resulting magnetospheric structure will be complex, and 
cannot be described in detail by a global _Ra, as testified by 
the complex Ha variations (Shultz et al. in prep). 



A35 HD 136504 (eLup) 

e Lup is an eccentric double-lined spectroscopic binary (P — 
4.6 d, e = 0.28; Uytterhoeven et al. 2005), with similar com- 
ponents. The magnetic field measurements found in the lit- 
erature (Hubrig et al. 2009; Shultz et al. 2012) do not specify 
which component is magnetic, but follow-up MiMeS obser- 
vations show that the field is associated with the primary 
star and that the published longitudinal field measurements 
are underestimated by a factor of two. 

Uytterhoeven et al. (2005) also found a possible period 
of 1.2 d for the primary. This would lead to W = 0.38 and 
R K = 1.9R,. 



A40 HD 186205 

The SIMBAD database gives a spectral type of B5, although 
Walborn (1975) classified it as B2Vp He-strong. The ef- 
fective temperature determinations are varied: 17 kK (Zbo- 
ril & North 2000) and 23.5 kK (Lee & Obrien 1977). An 
analysis of a MiMeS observation suggests T e s ~ 20 kK and 
logg = 4.0. 



A42 HD 200775 (V3780Cep) 

HD 200775 is a double-lined spectroscopic binary with a pe- 
riod ~ 4yr (e = 0.3) with components of similar tempera- 
tures. The magnetic field is associated with the sharp-lined 
primary star (Alecian et al. 2008a). 



A48 HD 58260 

According to the various parameter determinations reported 
by Cidale et al. (2007) and Bohlender (1989), we opt for 
T cff = 20 kK and log 5 = 3.55. 

Pedersen (1979) report a possible period of 1.657 d, 
based on spectrophotometry of the Hel A4026 A line, with 
a variation of the order of 0.01 mag. Bychkov et al. (2005) 
phased the dozen available longitudinal field measurements 
with this period, however the amplitude of the variation is 
relatively small. We therefore use the vsini as the lower 
limit on the equatorial velocity. The 1.6-d period, if con- 
firmed, would yield W = 0.55 and Rk ~ 1.5R,, and an 
inclination angle of the rotation axis near zero. 



A49 HD 36485 (5 Ori C) 

S Ori C is a double-lined spectroscopic binary with a 30-d 
period (e = 0.32). Leone et al. (2010) determined T e g — 
20 kK and M = 7Mq for the magnetic primary, and T e & = 
10 kK and M = 2.8 Mq for the secondary star (with a AV 
difference of 1.3 mag). 

Leone et al. (2010) determined a dipolar field strength 
between 7.3 and 12 kG. We use a mean value of 10 kG for 
our calculations. 



A50 HD 208057 (16 Peg) 

The star was reported once to display Ha emission (Merrill 
& Burwell 1943), but there exists no confirmation of this 
emission (Henrichs et al. 2009). 



A52 HD 25558 (40Tau) 

A magnetic field detection was reported by Hubrig et al. 
(2009), but was refuted by Bagnulo et al. (2012) based on a 
re-analysis of the same dataset. However, a weaker field was 
detected with MiMeS observations. 



A53 HD 35298 

Although Bychkov et al. (2005) reported a longitudinal field 
with extremum at 3kG, Yakunin et al. (2011) found some 
larger values up to 5kG but also reported large variation 
in the field determined with lines from different elements. 
The photometric period used by Bychkov et al. (2005) was 
confirmed by additional longitudinal field measurements at 
the Dominion Astrophysical Observatory, with an extremum 
around 5kG (D. Bohlender, priv com). We opt for a conser- 
vative lower limit on the dipolar strength of 9kG. 



A55 HD 142990 (V 913 Sco) 

There are many different effective temperature determina- 
tions in the literature, ranging from 16.5 to 18.5 kK (e.g. 
Cidale et al. 2007). We use 17.5 kK as a mean value and 
assume \og(g) = 4.0. 



A56 HD 37058 (V 359 Ori) 

Glagolevskij et al. (2007) determined T e g = 17 kK and 
log(g) = 3.80 from spectral fitting. Landstreet et al. (2007) 
found a higher temperature of 20 kK. We however prefer the 
lower temperature given that the star is He-weak. 

A period of 14 d was reported by Pedersen (1979) based 
on spectrophotometry of He lines. On the other hand, By- 
chkov et al. (2005) phased the sparse longitudinal field mea- 
surements from the literature with a period of 1.022 d. An 
additional measurement was taken by Bagnulo et al. (2006), 
but the period is not precise enough to test the phasing. 
We therefore use the conservative 14 d period. The shorter 
period would lead to W = 0.57 and Rk = 1.5i?». 
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A57 HD 35502 

HD 35502 is a hierarchical spectroscopic triple system with 
a broad-lined magnetic B-type primary (HD 35502 A) and a 
companion (HD 35502 Bab) composed of two sharp-lined A- 
type stars (Bohlender et al. in prep). Borra (1981) lists a pos- 
sible period of 1.7 d, however this period is not compatible 
with the new longitudinal field measurements (P = 0.85 d). 

A60 HD 61556 (HR2949) 

HD 61555/6 is a visual pair, the light of which is combined 
in the Hipparcos identifier HIP 37229 (Rivinius et al. 2003). 
Both components were observed separately in the context of 
the MiMeS Project, and a magnetic field was detected only 
for HD 61556. 

A61 HD 175362 (Wolff's star) 

There is a large scatter of effective temperature determina- 
tions in the literature, from 14 to 17 kK. We adopt the tem- 
perature (T e ff = 15 kK) derived by Leone & Manfre (1997) 
and logg = 4.0. 

A62 HD 105382 (HR 4618) 

Although often classified at a Be star, Briquet et al. (2001) 
have shown that this is not the case and the Be classifica- 
tions probably resulted from accidentally observing the very 
nearby, well known Be star <5Cen. 



